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1.  Introduction.  Spaceborne  imaging  of  the  Earth’s  surface  by  synthetic  aperture  radars  (SAR)  may 
be  adversely  affected  by  the  ionosphere  that  causes  distortions  of  the  signals  emitted  and  received  by  the 
radar  antenna.  In  Section  2  of  the  report,  we  analyzed  those  distortions  for  the  inhomogeneous  ionosphere 
described  by  the  cold  plasma  model.  Based  on  the  analysis,  we  can  conclude  that  the  deterioration  of  SAR 
images  is  due  to  the  mismatch  between  certain  parameters  of  the  actual  received  signal,  which  is  slowed  down 
by  the  temporal  dispersion  in  the  ionosphere,  and  the  corresponding  parameters  of  the  matched  filter,  which 
is  taken  as  if  the  propagation  was  unobstructed.  Consequently,  to  improve  the  quality  of  the  images  the  filter 
must  be  corrected.  However,  to  get  the  appropriate  correction  one  needs  to  know  some  key  characteristics  of 
the  ionosphere  precisely  at  the  time  and  place  the  image  is  taken.  To  obtain  those  characteristics,  we  propose 
to  probe  the  terrain,  and  hence  the  ionosphere,  on  two  distinct  carrier  frequencies.  We  also  investigate  the 
performance  of  the  matched  filters  that  were  corrected  this  way,  and  show  that  the  final  quality  of  the 
images  indeed  improves.  The  exposition  of  the  material  in  Section  2  follows  our  recent  paper  [ST10],  see 
also  [Tsy09a]. 

The  start-stop  approximation  is  a  standard  tool  for  processing  radar  data  in  synthetic  aperture  imaging. 
It  assumes  that  the  antenna  is  motionless  when  a  pulse  is  emitted  and  the  scattered  signal  received,  after 
which  the  antenna  moves  to  its  next  sending/receiving  position  along  the  flight  track.  However,  when  the 
antenna  is  mounted  on  a  satellite,  as  opposed  to  an  airplane,  its  relatively  high  speed  raises  at  least  two 
questions.  The  first  one  is  whether  the  image  may  be  affected  by  the  actual  displacement  of  the  antenna 
during  the  pulse  round-trip  time  between  the  orbit  and  the  Earth’s  surface.  This  displacement,  in  fact,  can 
be  rather  large.  Nonetheless,  by  analyzing  the  corresponding  generalized  ambiguity  function  of  a  SAR  sensor 
(see  Section  3  of  the  report)  we  show  that  in  practice  this  issue  can  be  disregarded.  The  second  question  is 
related  to  the  Doppler  frequency  shift,  which,  again,  is  larger  for  spaceborne  radars  than  for  airborne  radars. 
In  the  early  SAR  studies,  this  frequency  shift  provided  a  venue  for  understanding  the  azimuthal  resolution 
of  a  radar.  However,  in  a  more  rigorous  analysis  based  on  the  generalized  ambiguity  function,  the  Doppler 
effect  is  typically  left  out  of  consideration.  In  Section  3  of  the  report,  we  also  show  that  for  the  image  to 
stay  largely  unaffected  by  Doppler,  the  frequency  shift  must  be  included  in  the  definition  of  a  matched  filter. 
Otherwise,  there  will  be  a  geometric  shift  (translation)  of  the  entire  imaged  scene  from  its  true  position,  and 
there  may  also  be  a  slight  deterioration  of  the  image  sharpness  (contrast).  The  exposition  of  the  material  in 
Section  3  follows  our  paper  [Tsy09b]. 

2.  Dual  carrier  probing  for  spaceborne  SAR  imaging. 

2.1.  Background.  Synthetic  aperture  radar  (SAR)  is  a  coherent  imaging  technology  which  is  based 
predominantly  on  the  phase  information  of  the  emitted  and  received  waves  (as  opposed,  say,  to  conventional 
photography,  which  is  based  on  the  amplitude).  In  particular,  spaceborne  SAR  imaging  of  the  Earth’s 
surface  involves  sending  specially  shaped  pulses  of  radio  waves  from  the  satellite  to  the  ground,  receiving  the 
scattered  response,  and  processing  it  by  means  of  a  matched  filter  [CheOl,  CB09].  In  doing  so,  the  received 
signal  can,  in  some  approximate  sense,  be  interpreted  as  a  Fourier  transform  of  the  ground  reflectivity 
function,  whereas  filtering  that  yields  the  image  is  analogous  to  the  inverse  transform.  The  filter  is  called 
matched  because  its  phase  is  supposed  to  match  that  of  the  received  signal. 

When  radar  signals  travel  back  and  forth  from  the  satellite  to  the  surface,  they  propagate  through  the 
Earth’s  ionosphere  and  hence  undergo  temporal  dispersion.  This  causes  delays  in  arrival  times.  If  the  filter 
does  not  take  those  delays  into  account,  a  mismatch  of  the  phase  occurs  and  accordingly,  the  quality  of  the 
image  deteriorates  [Tsy09a,ST10].  The  extent  of  deterioration  depends  on  the  carrier  frequency  of  the  radar, 
or  more  precisely,  on  its  ratio  to  the  Langmuir  frequency  of  the  ionospheric  plasma.  The  higher  the  carrier 
frequency  the  less  prone  the  image  to  ionospheric  distortions.  Therefore,  many  modern  spaceborne  SAR 
sensors  operate  in  the  microwave  band  on  the  frequencies  of  around  10  GHz  or  even  higher.  For  example,  the 
new  European  satellite  TerraS  AR-X  operates  in  the  X-band,  on  the  frequencies  between  7  GHz  and  11.2  GHz. 

However,  there  is  also  substantial  interest  toward  spaceborne  SAR  technology  on  lower  carrier  frequen¬ 
cies,  in  the  VHF  or  UHF  band,  i.e. ,  in  the  range  of  hundreds  of  megahertz.  Radars  of  this  type  have  a  better 
surface  penetrating  capability.  On  the  other  hand,  their  images  are  much  more  sensitive  to  ionospheric 
distortions.  Consequently,  special  means  need  to  be  developed  to  compensate  for  those  distortions. 

In  [Tsy09a,ST10],  we  have  analyzed  the  distortions  for  the  case  of  SAR  imaging  by  means  of  a  scalar  field 
that  propagates  through  the  inhomogeneous  ionosphere  described  by  the  cold  plasma  model.  We  have  also 
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proposed  to  probe  the  terrain,  and  hence  the  ionosphere,  on  two  distinct  carrier  frequencies.  This  allows  one 
to  determine  some  key  ionospheric  parameters  that  control  the  propagation,  and  then  correct  the  matched 
filter  accordingly.  We  have  also  investigated  the  performance  of  the  matched  filters  with  the  corrections 
introduced  using  dual  carrier  probing.  Our  main  conclusion  is  that  the  new  filters  are  capable  of  completely 
eliminating  the  distortions  of  SAR  images  due  to  the  deterministic  part  of  the  charged  particles  content  in 
the  ionosphere,  including  the  inhomogeneous  case.  Removing  the  distortions  due  to  the  random  part  of  the 
charged  particles  content  will  be  subject  of  the  future  work. 

We  emphasize  that  the  corrections  shall  be  obtained  precisely  for  the  time  and  place  the  image  is  taken. 
The  reason  is  that  the  Earth’s  ionosphere  is  a  very  “lively”  medium;  its  characteristics  depend  on  many 
factors,  such  as  time  of  the  day,  time  of  the  year,  geographic  location,  level  of  Solar  activity,  etc.,  and  may 
change  rapidly,  e.g.,  considerably  faster  than  the  period  of  one  revolution  of  the  satellite  around  the  Earth. 

This  part  of  the  report  is  organized  as  follows.  In  Section  2.2,  we  describe  the  physical  model  and 
introduce  the  key  assumptions.  In  Section  2.3,  we  build  the  generalized  ambiguity  function  of  a  SAR,  which 
provides  a  key  tool  for  analyzing  the  radar  resolution.  In  Section  2.4,  we  discuss  the  performance  of  a 
spaceborne  SAR  sensor  in  the  case  when  no  attempt  is  made  to  mitigate  the  distortions;  this  creates  a 
foundation  for  the  subsequent  discussion  of  a  corrected  matched  filter.  In  Section  2.5,  we  introduce  the  dual 
carrier  probing  as  a  venue  for  obtaining  the  parameters  of  the  ionosphere  needed  for  correcting  the  filter. 
In  Section  2.6,  we  discuss  the  performance  of  the  radar  once  the  corrected  filter  has  been  implemented, 
and  show  that  the  resolution  and  sharpness  of  the  images  indeed  improve.  Finally,  Section  2.7  presents  a 
summary  of  what  has  been  accomplished,  and  identifies  directions  for  future  work. 

2.2.  The  model.  For  the  following  discussion  we  adopt  the  same  model  as  in  [Tsy09a,  ST10].  We 
analyze  the  standard  (as  opposed  to  interferometric)  SAR  images  [FL99].  We  assume  that  the  radar  does 
not  measure  the  polarization,  and  hence  consider  the  propagation  of  scalar  quantities  only.  All  the  tar¬ 
gets  are  assumed  deterministic  and  dispersionless;  scattering  of  radar  signals  off  the  targets  is  linearized 
and  interpreted  using  the  first  Born  approximation.  The  start-stop  approximation  is  employed  to  describe 
the  synthetic  antenna;  its  applicability  to  spaceborne  SAR  sensors  (as  opposed  to  airborne)  is  justified 
in  [Tsy09b],  see  also  Section  3  of  this  report. 

The  ionosphere  is  modeled  as  a  layer  of  inhomogeneous  dilute  plasma,  for  which  the  mean  concentration 
(number  density)  of  electrons  depends  on  the  altitude  but  does  not  depend  on  the  horizontal  coordinates. 
The  latter  assumption  holds  provided  that  the  characteristic  distances  along  the  Earth’s  surface  are  not 
very  large.  We  disregard  the  effect  of  the  magnetic  field  of  the  Earth  on  mean  electron  concentration  and 
assume  that  the  ionosphere  is  isotropic.  This  implies,  in  particular,  that  the  Faraday  rotation  is  not  taken 
into  account;  of  course,  for  scalar  propagating  fields  (no  polarization)  it  is  not  considered  anyway. 

On  top  of  its  mean  value,  the  electron  number  density  has  a  stochastic  component  due  to  the  ionospheric 
turbulence,  which  is  a  quasi-homogeneous  random  field  [MY71,RKT89a].  The  effect  of  the  Earth’s  magnetic 
field  on  the  ionospheric  turbulence  is  disregarded,  and  the  turbulence  is  assumed  isotropic. 

The  propagation  of  electromagnetic  waves  in  the  ionosphere  is  analyzed  with  the  help  of  the  cold  plasma 
approximation  [Gin64],  which  requires  that  the  phase  speed  of  the  waves  be  much  faster  than  the  thermal 
speed  of  the  electrons.  The  cold  plasma  approximation  is  equivalent  to  disregarding  spatial  dispersion  and 
taking  into  account  only  temporal  dispersion  of  radio  waves  in  the  plasma  [MM91].  The  ionospheric  plasma 
is  also  assumed  lossless,  i.e. ,  its  Ohm  conductivity  is  disregarded.  This  means  that  the  effective  collision 
frequencies  between  different  species  of  particles  in  the  plasma  are  much  lower  than  the  Langmuir  frequency. 
Moreover,  the  propagation  is  assumed  linear,  because  the  power  of  signals  emitted  by  spaceborne  antennas 
is  typically  much  lower  than  that  needed  for  setting  off  the  nonlinear  effects,  see  [GG60,  Gur78,  Gur07l,  as 
well  as  [Gin64,  Ch.  VIII]. 

We  use  the  classical  Helmholtz  theorem  [MF53]  to  partition  the  overall  electric  field  into  the  longitudinal 
and  transverse  components.  Then,  from  the  full  Maxwell’s  equations  and  the  equation  of  motion  of  electrons 
in  the  electric  field  one  can  derive  the  Klein-Gordon  equation  (see  [Tsy09a,  App.  A])  that  governs  the 
propagation  of  high  frequency  transverse1  electromagnetic  waves  in  dilute  plasma: 


d2E± 

dt2 


c2  AE±  +  oj2cE±  —  0. 


(2.1) 


1  The  corresponding  longitudinal  oscillations  are  known  as  Langmuir  waves,  they  have  zero  group  velocity  [MM91]. 
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The  quantity  wpe  in  equation  (2.1)  is  the  Langmuir  frequency,  or  plasma  electron  frequency.  It  characterizes 
temporal  responses  of  the  plasma  and  is  given  by: 


— 


1 ATre2Ne 


me 


(2.2) 


where  e  and  me  are  the  charge  and  mass  of  the  electron,  and  Ne  is  the  electron  number  density.  The 
behavior  of  wpe  in  the  ionosphere  follows  that  of  Ne;  it  has  an  altitude-dependent  deterministic  component 
and  a  stochastic  component  due  to  the  turbulence.  Typical  values  of  the  Langmuir  frequency  in  the  Earth’s 
ionosphere  range  between  3 MHz  and  15 MHz. 

Unlike  in  a  vacuum,  the  propagation  speed  of  electromagnetic  waves  in  the  ionosphere  depends  on 
the  frequency.  This  phenomenon  is  known  as  temporal  dispersion;  it  is  due  to  the  presence  of  the  last, 
non-differentiated,  term  in  equation  (2.1).  The  dispersion  relation  for  the  Klein-Gordon  equation  reads: 

w2=wpe+c2fc2.  (2.3) 

Dispersion  (2.3)  is  anomalous,  i.e. ,  it  is  the  short  waves  that  are  weakly  dispersive,  whereas  the  long  waves 
are  subject  to  stronger  dispersion.  Indeed,  the  phase  and  group  velocity  of  the  propagation  are  given  by: 

vph=c(l  +w2e/c2fc2)5  (2.4) 

and 

ugr  =c(l  Tw2e/c2fc2)-3  ,  (2.5) 

so  that  the  shorter  the  wave,  (the  larger  the  k)  the  closer  both  velocities  to  the  non-dispersive  limit  v  =  c. 
In  the  Cartesian  coordinates  we  can  consider  equation  (2.1)  for  individual  held  components: 

—  ^c2AE  +  uJ2peE  =  0.  (2.6) 

Hereinafter,  we  will  use  the  scalar  governing  equation  (2.6)  to  describe  the  propagation  of  electromagnetic 
waves  in  the  ionosphere.  For  a  comprehensive  account  of  the  propagation  of  radio  waves  in  plasma  we  refer 
the  reader  to  the  monographs  [Gin64]  and  [Bud85]. 

To  characterize  the  ionospheric  turbulence,  we  write  the  electron  number  density  as  follows: 


Ne  =  (Ne)  +  fi(x), 


(2.7) 


where  the  angular  brackets  ( • )  denote  the  expected  value  (mean)  and  \j  represents  the  fluctuations:  (/z)  =  0. 
In  the  simplest  case  of  constant  (Ne),  n(x)  is  a  homogeneous  and  isotropic  random  held,  and  its  correlation 
function  depends  only  on  the  distance  r  =  \x±  —  a^|  and  not  on  the  individual  locations  x1  and  x2: 

V(xltx2)  =  (fi(x1)/u(x2))  =  V(\x!  —  x2\)  =  V(r).  (2.8) 


In  reality,  however,  (Ne)  is  a  function  of  the  altitude  h.  At  the  same  time,  the  ratio 


M  = 


vM 

(Ne) 


(2.9) 


is  assumed  altitude  independent.  A  typical  numerical  value  of  M  is  5-  1CU3,  and  in  extreme  situations  it  may 
reach  10-1  [Arm05].  Hence,  \J (fi2)  also  depends  on  h,  which  makes  )  a  quasi-homogeneous  (rather  than 
truly  homogeneous)  random  held.  While  in  this  case  we  still  keep  the  same  notation  (2.8)  for  the  correlation 
function,  in  fact  we  have  V  =  V(h,r),  where  the  dependence  on  h  (altitude)  is  slow  and  the  dependence  on 
r  (local  variable)  is  fast.  Hereafter,  we  will  use  the  exponential  correlation  function 

V(r)  =  Ctt2  qpe~q°r . 


(2.10) 
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It  accounts  only  for  the  short  range  phenomena  in  the  ionosphere  and  has  the  spectrum  [MY75,  Sec.  12.1] 


V(q) 


c 

(i  +  q2  /  <?o)K  5 


where  k  =  2  and  C  =  const. 


(2.11) 


In  formulae  (2.10)  and  (2.11),  we  denote  =  r q  ,  where  r0  is  the  correlation  length  or  outer  scale  of 
turbulence.  The  latter  may  range  between  1  km  and  10 km  according  to  different  sources  in  the  literature. 
An  advantage  of  using  the  correlation  function  (2.10)  is  that  one  can  easily  make  sure  that  it  satisfies  a 
necessary  and  sufficient  condition  of  ergodicity  given  in  [MY71,  Sec.  4.7]  (a  1938  theorem  by  E.  E.  Slutskii): 


s 

&  I  /  v{s)d‘ = °- 

0 


(2.12) 


Ergodicity  implies  that  statistical  means  can  be  substituted  for  spatial  averages  of  /z  whenever  averaging  is 
done  over  sufficiently  long  intervals.  In  that  regard  we  note  that  if  a  random  field  is  ergodic,  then  it  must 
have  (/j,)  =  const,  because  the  spatial  average  of  /j  does  not  depend  on  the  location,  whereas  its  statistical 
mean  (fi)  may,  generally  speaking,  be  a  function  of  the  coordinates  [MY71,  Sec.  4.6].  In  our  particular  case, 
however,  we  have  (/z)  =  0,  which  means  that  ergodicity  can  be  employed  even  when  the  random  field  /i  is 
quasi-homogeneous  rather  than  truly  homogeneous  (because  its  variance  depends  on  the  altitude). 

It  should  also  be  mentioned  that  for  the  true  Kolmogorov- type  turbulence  we  have  k  =  11/6  rather  than 
k  =  2  in  formula  (2.11),  see  [RKT89a,  Ch.  I].  From  the  standpoint  of  practice,  the  difference  between  the 
two  cases  is  that  n  =  11/6  allows  one  to  take  into  account  some  long  range  phenomena  in  the  ionosphere. 
As  we  have  indicated  in  [Tsy09a],  this  is  a  separate  interesting  subject  for  the  future  study. 

For  the  purpose  of  obtaining  specific  quantitative  estimates,  we  will  use  the  same  typical  values  of  the 
parameters  as  we  have  exploited  in  [Tsy09a,ST10].  A  reference  value  of  the  mean  electron  concentration  will 
be  taken  as  (Ne)  =  10 6 cm-3  (in  fact,  (Ne)  =  ( Ne(h ))),  which  corresponds  to  the  F-layer  of  the  ionosphere 
and  yields  wpe  «  9 MHz.  The  radar  carrier  frequency  will  be  assumed  equal  to  1  GHz,  which  is  the  lower 
bound  for  the  microwave  band,  and  which  yields  the  wavelength  A  =  30cm.  According  to  [Bro93],  it  is 
very  difficult  to  achieve  high  resolution  of  spaceborne  S AR  images  for  carrier  frequencies  under  1  GHz  if  no 
attempt  is  made  to  reduce  the  distortions  due  to  the  ionosphere.  A  typical  one-way  propagation  distance 
between  the  satellite  and  the  ground  will  be  taken  as  Rq  ss  lOOOfcm,  and  the  orbit  altitude  H  may  be  about 
500 km.  or  higher.  Other  important  parameters  will  be  introduced  as  needed. 

2.3.  The  generalized  ambiguity  function.  The  generalized  ambiguity  function  offers  a  convenient 
way  of  assessing  the  radar  performance.  It  is  basically  the  image  of  a  point  target,  i.e.,  that  of  a  delta-type 
scatterer.  In  the  ideal  world,  this  image  would  be  a  delta-function  as  well.  In  reality,  however,  it  is  smeared 
out  (by  the  nature  of  the  radar  signal  processing,  even  when  there  is  no  ionosphere  and  no  dispersion),  and 
the  extent  of  this  smearing  provides  the  limits  of  radar  resolution. 

The  interrogating  pulses  emitted  by  the  antenna  of  a  SAR  are  taken  as  linear  upchirps2  of  the  form: 

P{t)  =  A(t.)eiuiot,  where  A{t)  =  Xr(t)eiat2 .  (2.13) 


In  formula  (2.13),  yT(t)  is  the  indicator  function  of  the  interval  of  duration  r: 


Xr(t) 


1,  t  £  [— t/2,  t/2], 
0,  otherwise, 


and  a  =  B/(2t),  where  B  is  the  bandwidth  of  the  chirp.  Accordingly,  the  instantaneous  frequency  of  the 
chirp  [Che01,CB09]  is  given  by 


Bt 

u{t)  =  w0  H - ;,  te  [-t/2,  t/2], 

T 


(2.14) 


2  High-range  resolution  waveforms,  see  [CB09]. 


AFOSR  Grant  No.  FA9550-07-1-0170.  Final  Performance  Report 


7 


where  ui o  is  the  center  carrier  frequency.  The  modulating  function  A(t)  in  formula  (2.13)  is  assumed  slowly 
varying  compared  to  the  fast  carrier  oscillation  elUot.  A  typical  duration  of  the  pulse  in  actual  SAR  systems 
may  be  r  ~  5  •  10_5sec,  a  typical  interval  between  two  consecutive  pulses  ~  5  •  10_4sec,  and  the  bandwidth 
in  formula  (2.14)  is  B  ~  10 MHz  (it  may  be  higher). 

Hereafter,  the  radar  is  assumed  to  operate  in  the  stripmap  mode,  see  [CGM95],  when  the  antenna  points 
in  a  fixed  direction  relative  to  that  of  the  satellite  motion.3  Hence,  the  footprint  of  the  beam  emitted  by  the 
antenna  sweeps  a  strip  on  the  Earth’s  surface  parallel  to  the  flight  track,  i.e. ,  to  the  orbit. 

To  obtain  the  image,  the  antenna  emits  a  series  of  pulses  (2.13)  when  moving  along  the  orbit.  Those 
pulses  travel  to  the  Earth’s  surface,  get  scattered  off  the  ground,  and  propagate  back  toward  the  satellite 
where  they  are  received  by  the  same  antenna  (in  the  case  of  a  monostatic  SAR  [CGM95]).  In  doing  so, 
clearly,  both  the  emission  and  the  reception  of  pulses  are  done  by  a  moving  antenna  (which  gives  rise  to  the 
Doppler  frequency  shift)  and  moreover,  the  satellite  travels  a  certain  distance  along  the  orbit  during  the  pulse 
round-trip  time  between  the  antenna  and  the  ground.  However,  in  a  simplified  framework  of  the  start-stop 
approximation  those  effects  are  disregarded.  Instead,  we  assume  that  the  pulse  is  emitted  and  the  scattered 
response  received  while  the  antenna  is  at  standstill  at  one  and  the  same  position,  after  which  it  moves  to 
its  next  sending/receiving  position  along  the  orbit.  This  assumption  has  been  used  extensively  in  the  SAR 
literature  as  it  renders  the  corresponding  analysis  much  easier.  It,  however,  needs  to  be  justified,  especially 
when  the  antenna  moves  fast,  such  as  when  it  is  mounted  on  a  satellite  as  opposed  to  an  airplane.  A  full 
justification  for  the  use  of  start-stop  approximation  for  spaceborne  SAR  imaging  can  be  found  in  [Tsy09b], 
see  also  Section  3  of  this  report. 


2.3.1.  Non-dispersive  propagation.  We  first  derive  the  generalized  ambiguity  function  for  the  case 
of  unobstructed  propagation  with  the  speed  c.  In  this  case,  the  field  is  governed  by  the  standard  d’Alembert 
equation,  which  is  obtained  from  the  Klein-Gordon  equation  (2.6)  by  dropping  the  last  term,  i.e.,  by  setting 
wpe  =  0.  Suppose  that  the  antenna  is  a  motionless  point  source  located  at  a;  €  K3.  Then,  the  propagating 
held  due  to  the  emitted  chirp  (2.13)  is  given  by  the  standard  retarded  potential  of  the  d’Alembert  operator: 


1  P(t  —  \z  —  x\/c) 
47t  |  z  —  x  | 


(2.15) 


Let  us  assume  that  the  imaged  terrain,  which  is  also  motionless,  is  characterized  by  the  variable  refraction 
index  n  =  n(z).  Under  the  first  Born  approximation  [BW99,  Sec.  13.1.4],  scattering  is  linearized  so  that  the 
terrain  is  interpreted  as  a  secondary  waves’  source  due  to  the  incident  held  <p(t,  z)  of  (2.15): 


1  -  n2(z)  d2ip  def  (  ^  d2(p 

^  ~  u[z)~di 2"' 

Consequently,  the  scattered  held  at  the  location  x  of  the  antenna  and  the  moment  of  time  t  is  given  by  the 
Kirchhoff  integral: 

=  i  HI  (2-16) 

\x— z\<ct 


As  the  amplitude  A(t)  in  (2.13)  is  slowly  varying,  it  can  be  left  out  when  differentiating  the  incident  held 
(2.15)  for  substitution  into  (2.16),  which  yields: 


d2<p  ^  ' ,2 

dt 2  ’  47T  I  z  —  x 


Wq  P{t~  \Z  -  x\/c) 


(2.17) 


Consequently, 


*!’{t,x)  «  - 


fff  "(z) 

16tt2  JJJ  \x  -  z | 

\x— z\<ct 

v(z) 


16tt2  JJJ  \x  -  z\ 

\x  —  z\<ct 


,P  (t  —  2\x  —  z\/c)  dz 


;A  (t  -  2\x  -  z\/c)  eiuo(t-2\x-z\/c)  dz 


(2.18) 


3 When  this  direction  is  normal  to  the  flight  track,  the  imaging  is  referred  to  as  broadside. 
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According  to  (2.18),  the  scattered  field  ip(t,  x )  can  be  interpreted  a  result  of  application  of  a  Fourier  integral 
operator  (FIO)  to  the  ground  reflectivity  function  i/(z),  see,  e.g.,  [CN04,  NC04,  CB08,  CB09].The  FIO  is 
approximately  inverted  by  applying  a  matched  filter  to  0  and  accumulating  the  information  due  to  multiple 
interrogating  pulses  (2.13)  emitted  from  and  received  by  the  antenna  at  different  locations  on  the  orbit.  This 
procedure  is  similar  to  application  of  the  adjoint  operator,  which  would  have  coincided  with  the  true  inverse 
if  the  mapping  (2.18)  was  a  standard  Fourier  transform. 

The  matched  filter  is  defined  as  follows.  Assume  that  there  is  a  point  scatterer  at  the  reference  location 
y,  then  the  resulting  field  at  (t,  x)  is  obtained  by  substituting  i/(z)  =  S(z  —  y)  into  formula  (2.18): 


0i  (t,x)  = 


Up  P(t-2\x-  y\/c) 
167T2  \x  —  y  |2 


(2.19) 


The  filter  is  essentially  a  complex  conjugate  of  0 1  given  by  (2.19);  for  simplicity,  the  constant  factor 
— Wq/167 t2,  as  well  the  entire  denominator,  which  is  a  slowly  varying  function  (compared  to  the  fast  os¬ 
cillation  elUJot),  are  disregarded.  What  remains  is  merely  P(t  —  2\x  —  y\/c),  where  the  overbar  denotes 
complex  conjugation.  The  idea  of  building  the  filter  this  way,  i.e.,  the  idea  of  matching,  is  to  have  the  large 
phase  cancel  in  the  exponent,  see  formula  (2.13).  The  application  of  this  filter  yields  a  single-look  image: 


I(y)=  j  P(t-2\x  -  y\/c)ip{t,x)dt 


(2.20) 


107T2 


P(t  —  2\x  —  y\/c)P(t  —  2\x  —  z\/c)dt  ■ 


v(z) 


X  -  z 


■  dz , 


\x—  z\<ct 


W{y,z) 


where  we  have  changed  the  order  of  integration  after  substituting  expression  (2.18)  for  ip(t,  x).  The  interior 
integral  W(y,z)  in  formula  (2.20)  is  called  the  point  spread  function,  see  [CheOl,  CB09].  Up  to  a  slowly 
varying  denominator,  the  point  spread  function  W(y,z)  yields  a  single-look  image  of  the  point  scatterer 
located  at  z,  i.e.,  it  is  the  field  due  to  a  unit  magnitude  delta- function  at  z  processed  with  the  matched  filter 

As,  however,  has  been  mentioned,  the  full  SAR  image  is  built  using  multiple  looks,  when  the  antenna 
emits  a  sequence  of  chirps  (2.13)  from  different  locations  on  the  orbit.  It  is  precisely  this  approach  based 
on  incorporating  the  information  from  multiple  looks  that  enables  drastic  improvements  in  the  quality  of 
the  image  compared  to  the  single  look  strategy  [CheOl,  CB09].  Let  us  therefore  consider  a  sequence  of 
emitting/receiving  times  and  locations  ( tn,xn ).  For  each  n,  we  build  the  point  spread  function  following 
(2.20): 


Wn{y,z)  =  J  P(t  -  tn  -  2 |®n  -  y\/c)P(t  -  tn  -  2\xn  -  z\/c)dt.  (2.21) 

The  generalized  ambiguity  function  of  a  SAR  system  takes  into  account  the  information  from  multiple 
interrogating  pulses  by  summing  up  the  corresponding  contributions  (2.21): 

W(y,z)  =Y2^(z,xn)Wn(y,z).  (2.22) 

n 

The  factor  i?(z,  xn )  under  the  sum  in  (2.22)  determines  the  range  of  summation.  It  comes  from  the  directivity 
pattern  of  the  antenna,  because  the  actual  SAR  antenna  is  never  a  point  monopole.  The  quantity  i9(z,  xn) 
can  be  approximated  as  follows: 


d(z,xn) 


1,  if  the  target  z  is  in  the  beam  emitted  from  xn, 
0,  otherwise. 


In  [Tsy09a,  App.  B]  (see  also  [CheOl,  CB09])  we  have  shown  that  the  antenna  of  longitudinal  dimension 
L  emits  a  beam  of  angular  width  2X/L  provided  that  the  carrier  wavelength  A  =  27tc/wo  is  much  shorter 
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than  L.  Consequently,  the  longitudinal  size  of  the  antenna  footprint  on  the  ground  in  the  case  of  broadside 
imaging  is 


Ro  tan  «  R0^-, 

where  Ro  is  the  distance  from  the  location  of  the  antenna  xn  to  the  center  of  the  footprint.  Therefore, 


i?(z,  xn)  = 


1,  if  zi  -  R0j  <  x ?  <  zi  +  Raj 


L  ’ 


0,  if  x i  <  z\  -  Rqt  or  xi  >  zi  +  ^oji 


(2.23) 


where  the  subscript  “1”  denotes  the  coordinate  direction  along  the  flight  track  (orbit). 

The  range  of  summation  in  (2.22)  is  therefore  defined  by  taking  into  account  only  those  xn,  for  which 
d(z,  xn )  =  1,  see  (2.23).  Let  Aaq  be  the  distance  along  the  orbit  between  the  successive  emissions  of  pulses. 
We  can  also  take  zi  =  0  with  no  loss  of  generality.  Then,  the  inequality  <  a;”  <  translates  into 


—  R  <n<  R,  where  N  = 


2A  Rq 


Axi  L 


and  [  •  ]  is  the  integer  part.  Consequently,  we  recast  (2.21),  (2.22)  as 


N/2 

W(y,  z)  =  £  A(t-tn-2\y-xr'\/c)e2i“°\y-*n^ 

—n/2  J  (2.24) 

xA(t  -tn-  2|z  -  xn\ /c)e-2iu,^z~x^/cdt. 


Next,  we  change  the  integration  variable  from  t  to  t  —  tn  in  each  term  of  the  sum  (2.24)  and  realize  that 
neither  A(t  —  2| y  —  xn\/c)  nor  A(t  —  2|z  —  xn\ /c)  depends  on  n  explicitly,  except  for  the  dependence  via  xn. 
The  latter  is  weak,  because  for  large  Ro  both  \y  —  xn\  and  |z  —  xn\  are  slowly  varying  functions  of  xn,  and 
they  appear  as  arguments  of  another  slowly  varying  function,  A.  Hence,  the  factors  A  and  A  can  be  taken 
out  of  the  sum  (2.24),  and  xn  inside  A{  ■  )  and  A{  ■  )  can  be  replaced  by  x°  for  definiteness,  which  yields: 


f _  N/2 

W(y,z)  «  (  /  A{t-2\y  -  x°\/c)A(t  -  2|z  -  xa\/c)dt^j  ( 

y  AT  /O 


o2iui0(\y-xn\/ c-\z-xn\/ c) 


-N/2 


(2.25) 


WR(y,z) 


WA(y,z) 


Hence,  the  generalized  ambiguity  function  gets  approximately  split  into  the  product  of  the  range  factor: 

WR(y,z)  =  J  A{t-2\y  -  x°\/c)A{t-2\z  -  xa\/c)dt  (2.26) 

and  the  azimuthal  factor: 

N/2 

WA{y,z)  =  E  e2iu0(\y- Xn\fc- \z- xn|/c)_  (2.27) 

n—— N/2 

that  control  the  range  and  azimuthal  resolution  of  the  radar,4  respectively,  see  [ClieOl,  CB09],  in  the  case 
of  unobstructed  propagation  between  the  orbit  and  the  ground.  Indeed,  the  generalized  ambiguity  function 
(2.25)  is  the  image  of  a  delta-type  scatterer.  Hence,  if  we  could  make  it  equal  to  a  delta-function  as  well, 
W(y,z)  =  S(y  —  z),  then  the  radar  would  have  had  perfect  resolution.  In  reality  this  is  never  achieved, 
and  instead  W(y,  z)  has  a  peak  of  finite  height  and  finite  width  when  the  reference  point  y  approaches  the 
target  z.  The  sharper  (i.e. ,  the  narrower)  this  peak  the  better  the  resolution,  because  two  sharper  peaks  can 
be  told  apart  if  they  are  closer. 


4Resolution  is  a  capability  of  the  radar  to  distinguish  between  two  different  targets  located  a  certain  distance  apart. 
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2.3.2.  Dispersive  propagation  in  the  homogeneous  medium.  In  the  case  of  propagation  through 
the  ionosphere,  both  the  shape  and  the  speed  of  the  pulse  change  because  of  the  dispersion.  The  simplest 
case  to  analyze  is  that  of  a  homogeneous  ionosphere;  a  detailed  analysis  (based  on  the  Fourier  transform  and 
linearization)  is  provided  in  Appendix  A.  The  form  of  the  propagating  pulse  emitted  by  the  antenna  at  x  is 
given  by  equation  (A. 16)  [cf.  formula  (2.15)]: 


in{t  z)  =  ^  X\lV&(U°))ciun(t-\z-x\/v^(u0')') 

47r\z  —  X  | 


(2.28) 


where  A' (t)  =  Xr'(f)eiQ  4  , 

and  t'  and  a!  are  the  new  chirp  duration  and  rate  given  by  equations  (A.  12)  and  (A.  15),  respectively: 


/  ,  x  |  r>  I Z  X I  Wpe  B 

T  =  T  +  OT  =  T  +  2 - ^ - , 

C  Wg  Wo 


and 


/  ,  r  \z  X\  Wpe  B2 

a  =  a  +  da  =  a - 


St  B 


Wn  T2  a  2  T2 


B 

27 


1-^ 

T 


(2.29) 


(2.30) 


Let  us  emphasize  that  unlike  in  a  vacuum,  the  shape  of  the  pulse  keeps  changing  as  it  propagates  through 
the  plasma;  according  to  (2.29)  and  (2.30),  the  longer  the  distance  | z  —  x\  that  the  pulse  travels,  the  more 
it  gets  dilated  (r'  becomes  longer)  and  the  more  its  rate  decreases  ( a '  becomes  smaller). 

We  also  emphasize  that  as  formula  (2.28)  indicates,  the  chirp  itself,  i.e.,  its  slowly  varying  envelope  A' 
which  presents  a  physical  observable,  travels  with  the  group  velocity  ugr(wo)  that  corresponds  to  the  center 
carrier  frequency  u>o-  For  the  high  frequency  case  that  we  are  interested  in,  wo  wpe,  the  group  velocity 
(2.5)  can  be  linearized: 


^gr 


=  Ugr(w)  S3  C 


|  _  1  Wpe 


2  c2k2 


=  c  (l-l 


pe 


2  w2  —  w2e 


(2.31) 


At  the  same  time,  the  center  carrier  oscillation  per  se,  i.e.,  the  wave  elLUot,  travels  with  the  corresponding 
phase  velocity  uph(wo)  that  can  also  be  linearized: 


Vph  =  Uph(w)  S3  C 


1  Woe 

1  H - ^ 

2  c2k2 


=  C  !+2W^V 

V  2  w2  -  wpe 


1  + 


2  w2 


(2.32) 


Using  the  notions  of  the  group  and  phase  velocity  (2.31)  and  (2.32),  we  can  introduce  the  group  and  phase 
travel  times,  respectively,  in  the  homogeneous  ionosphere: 


Tgr(m,  z,  w)  — 

and  recast  formula  (2.28)  as  follows: 


ugr(w) 


and  Tph(x,  z,  w)  = 


Wh(w)  ’ 


m(t  z)  =  — T&r(X'  ' 

47t|z  —  x\ 


(2.33) 


(2.34) 


It  is  important  to  point  out  that  the  group  velocity  ugr  of  (2.31)  is  slower  than  the  speed  of  light  c,  as  it  is 
supposed  to  be  for  any  velocity  associated  with  the  propagation  of  physical  observables.  However,  the  phase 
velocity  uph  (2.32)  is  faster  than  the  speed  of  light  c,  because  the  dispersion  in  the  ionosphere  is  anomalous. 
Hence,  if  we  introduce  the  neutral  (non-dispersive)  travel  time  T  =  T{x,z)  =  \x  —  z |/c,  then,  according  to 
(2.31)  and  (2.32),  for  the  group  and  phase  travel  time  we  can  write: 


Tgr(x,  z,  w)  =  T(x,  z)  +  AT  and  Tph(x,  z,  w)  =  T(x,  z)  —  AT, 

where  AT  =  AT(x,  z,  w)  =  — - — 

c  2  w2 


(2.35) 
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In  other  words,  the  group  delay  and  the  phase  advance  in  formula  (2.34)  are  equal  to  one  another.  Given 
(2.35),  we  can  recast  formula  (2.29)  for  the  pulse  dilation  as  follows: 

TD 

6t  =  4AT(x,  z,  wo)  — .  (2.36) 

w0 


It  shall  also  be  noted  that  in  our  work  [Tsy09a],  we  used  a  simpler  semi-qualitative  approach  to  analyze 
the  propagation  of  chirps  in  the  ionosphere.  The  result  was  the  same  as  (2.34)  with  one  important  difference 
though:  Instead  of  the  phase  advance  that  characterizes  the  propagation  of  the  carrier  frequency  in  formula 
(2.34),  in  [Tsy09a]  we  have  incorporated  a  phase  delay  equal  to  the  group  delay.  This  could  not  affect 
the  analysis  of  the  range  resolution  as  it  is  based  on  the  propagation  of  envelopes,  see  formula  (2.26).  On 
the  other  hand,  this  could  potentially  affect  our  analysis  of  the  azimuthal  resolution  in  [Tsy09a],  as  the 
azimuthal  factor  of  the  generalized  ambiguity  function  (2.27)  includes  the  phase  information.  However,  as 
shown  in  [ST10],  our  previous  calculations  of  the  azimuthal  resolution  are  not  affected  either,  because  in  the 
linearized  high  frequency  framework  that  we  employ,  the  group  delay  and  the  phase  advance  are  equal  in 
magnitude,  see  formulae  (2.35). 

2.3.3.  Dispersive  propagation  in  the  inhomogeneous  medium.  In  the  actual  ionosphere,  the 
electron  number  density  Ne  is  not  constant,  and  this  is  going  to  affect  the  propagation  times  (2.33).  For  our 
subsequent  analysis,  we  will  assume  that  the  mean  electron  number  density  (Ne)  depends  on  the  altitude 
above  the  Earth’s  surface,  but  does  not  depend  on  the  horizontal  coordinates.  A  typical  dependence  of  the 
mean  electron  number  density  on  the  altitude  h  is  non- monotonic.  The  maximum  is  reached  in  the  F-layer 
somewhere  between  200fcm  and  300 km  above  the  Earth’s  surface,  and  the  characteristic  scale  h0  of  the 
variations  of  (Ne)  is  on  the  order  of  tens  of  kilometers,  see  [Gin64,  Ch.  VI].  Clearly,  ho  A,  where  A  is  the 
wavelength,  which  suggests  that  we  can  perhaps  use  the  approximation  of  geometrical  optics  to  analyze  the 
propagation  of  SAR  pulses  in  the  inhomogeneous  ionosphere. 

We  should  remember,  however,  that  the  total  electron  number  density  Ne  also  has  a  stochastic  component 
H(x),  see  formula  (2.7).  It  depends  on  all  spatial  coordinates,  and  can  be  interpreted  as  a  quasi-homogeneous 
random  field.  To  justify  the  use  of  geometrical  optics  for  the  study  of  pulse  propagation  through  the  turbulent 
ionosphere,  the  wavelength  A  must  be  much  shorter  than  the  characteristic  scale  of  turbulent  inhomogeneities. 
If  the  latter  is  taken  as  ro  (the  outer  scale  of  turbulence),  the  constraint  A  <C  ro  is  obviously  met.  There 
is,  however,  a  more  subtle  criterion  for  applicability  of  the  geometrical  optics.  The  characteristic  scale  of 
inhomogeneities  must  be  much  longer  than  the  size  of  the  first  Fresnel  zone  \/X Rz,  where  Rz  =  |z  —  x\  is 
the  propagation  distance  between  the  antenna  and  the  target.  The  quantity  \/\ Rz  comes  to  approximately 
540m  for  A  =  30cm  and  Rz  =  1000/cm,  which  is  roughly  2 Tq  according  to  [Arm05]  or  NrQ  according 
to  [BG88] .  Technically  speaking,  this  makes  the  geometrical  optics  a  borderline  approximation  for  the  class 
of  problems  we  are  considering.  It  is  known  however,  that  there  are  fewer  shorter  scale  inhomogeneities 
in  the  spectrum  of  ionospheric  turbulence  than  longer  scale  inhomogeneities,  which  still  leaves  the  main 
conclusions  of  geometrical  optics  valid  even  outside  its  formal  applicability  range,  see  [RKT89b,  Ch.  I], 

The  expressions  for  travel  times  in  the  inhomogeneous  ionosphere  that  replace  formulae  (2.33)  are  derived 
in  Appendix  B  for  the  deterministic  case  and  in  Appendix  C  for  the  stochastic  case.  The  final  expressions 
that  we  present  correspond  to  the  stochastic  formulae  (C.15)  and  (C.16);  they  reduce  to  the  deterministic 
formulae  (B.21)  and  (B.22)  if  /i  =  0: 


Tgv(x,  z,u>)  = 

Tph(x,z,u)  =  — 
c 

_  /  TT\ 

The  quantity  Ne  ;  in  formulae  (2.37)  is  defined  by  (C.8),  and  H  denotes  the  orbit  altitude.  Similarly  to 


L  +  - 


1  47re2  l  N, 


(H) 


2  ?new2  \  H 


L  - 


1  47re2  /  Nf 


(H) 


2  ?nPw2  \  H 


+  J  ^(s)ds 

\ 

n(s)d,s 


o 

Rz 


+ 


Rz 


(2.37a) 

(2.37b) 
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(2.35),  we  can  also  write: 


Tgr(x,  z,cj)  =  T(x,  z)  +  AT  and  Tph(®,  2,  w)  =  T(cc,  z)  —  AT, 


where  T(x,z)  =  — -  and  AT  =  AT(x,  z,  a;) 
c 


ftz  1  47re2  /  JVe(g) 
c  2  mew2  l  if 


1 


Rz 


y(s)ds 


M 


(2.38) 


Given  the  travel  times  (2.38),  we  can  evaluate  the  dilation  of  the  pulse  on  its  way  between  the  antenna  x 
and  the  target  z  using  formula  (2.36): 


ft  2ft  47TP2  ft 

5t  =  4AT(x,z,ujo)—  =  —  —  M.  (2.39) 

wo  c  meoj o  wo 

On  the  way  between  the  antenna  and  the  target  and  back,  the  dilation  (2.39)  doubles  because  the  travel 
distance  doubles.  The  pulse  rate  (2.30)  changes  accordingly. 

Finally,  we  can  modify  the  generalized  ambiguity  function  derived  in  Section  2.3.1,  see  formula  (2.25), 
by  substituting  the  actual  travel  times  (2.37)  instead  of  the  unobstructed  travel  time  | z  —  x°\/c  into  those 
factors  that  correspond  to  the  received  signal.  According  to  the  form  (2.34)  of  the  propagating  pulse,  the 
group  travel  time  Tgl-(x°,  z,  oj0)  of  (2.37a)  shall  be  substituted  into  WR{y,z)  and  the  phase  travel  time 
Tph(xn,  z,  wo)  of  (2.37b)  shall  be  substituted  into  WA(y,  z).  Then,  instead  of  formulae  (2.26)  and  (2.27)  we 
have: 


and 


W'R(y,z )  =  J  A(t-  2| y  -  x°\/c)A'2gT(t  -  2Tgr(x°,  z,  u0))dt 


(2.40) 


N/2 

W'A(y,z)=  ^  e2iw0(|y-®"|/c-Tph(a;'\z,w0))^ 

n=-AT/2 


(2.41) 


The  subscript  2 St  in  formula  (2.40)  indicates  that  the  round  trip  pulse  dilation  is  twice  the  St  of  (2.39). 

We  emphasize  that  in  formulae  (2.40)  and  (2.41)  the  factors  that  correspond  to  the  actual  field  received 
by  the  antenna  take  into  account  the  dispersion  of  radio  waves  in  the  ionosphere.  However,  the  factors  that 
represent  the  matched  filter  remain  the  same  as  in  the  non-dispersive  case.  This  creates  a  mismatch,  and 
as  has  been  shown  in  [Tsy09a],  it  is  precisely  this  mismatch  that  is  responsible  for  the  deterioration  of  the 
image.  In  the  next  section,  we  briefly  reproduce  the  corresponding  arguments  from  [Tsy09a].  It  will  also 
help  us  better  understand  the  mechanism  of  correcting  the  matched  filter. 

2.4.  Performance  of  the  radar  with  non-corrected  filter. 

2.4.1.  Range  resolution.  Integration  (2.40)  has  been  performed  in  [Tsy09a,  App.  C].  Up  to  a  factor 
of  magnitude  one,  integral  (2.40)  takes  the  form: 


W'R(y,z)  oc 


^min{r/2+2i?y/c,  t' /2-\-2TgI(x° 


ei{a' -a)t2  e4i(aRy  /c-a'  TgT(x°  ,z,u>0))t  ^ 


f  max{  —  t /2-\-2Ry  / c,  —  r'/2+2Tgr  (xQ  ,z 


(2.42) 


where  Ry  =  \y  —  x°\.  In  the  simplest  case,  we  disregard  the  quadratic  term  in  the  exponent  under  the  integral 
(2.42).  Then,  taking  the  limits  as  =Ft/2  +  2 Ry/c  and  changing  the  integration  variable:  t  =  u  +  2Ry/c,  we 
have: 


W'R(y,  z)  oc  t sine  (2 [a'Ry/c  —  a'TgI(x° ,  z,  w0)]t)  .  (2.43) 

def  * 

When  analyzing  radar  resolution,  the  sine  a;  =  function  is  what  one  typically  gets  instead  of  the 

“genuine”  (5-function,  see  [Che01,CB09].  The  sharper  the  central  peak  of  the  sine  (i.e. ,  its  main  lobe  between 
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the  first  zeros  on  either  side  of  the  maximum),  the  “closer”  it  is  to  the  5-function  and  hence  the  better  the 
resolution.  The  width  (or  sometimes,  semi-width)  of  the  main  lobe  of  the  sine  is  actually  taken  as  a  measure 
of  radar  resolution,  because  it  is  considered  the  distance  at  which  two  such  peaks  can  be  clearly  told  apart. 

In  the  non-dispersive  case  T(x°,  z,  w0)  =  Rz/c  and  a'  =  a.  Hence,  W'R{y ,  z)  reduces  to  Wr(p,  z)  and 
consequently,  the  sine  attains  its  maximum  when  Ry  =  Rz,  i.e. ,  when  the  distance  from  the  antenna  to 
the  reference  location  y  coincides  with  that  for  the  target  location  z.  The  resolution  is  then  controlled  by 
how  widely  the  sine  is  “spread  around”  its  central  maximum.  The  semi-width  of  the  main  lobe  of  the  sine 
corresponds  to  the  value  of  its  argument  equal  to  7 r  and  therefore,  for  the  full  width  we  find  A R  =  2-kc/ B. 
It  is  the  quantity  A R  that  is  interpreted  as  the  range  resolution,  i.e.,  the  minimum  distance,  at  which  the 
SAR  system  can  distinguish  between  two  point  targets  located  on  the  same  line  of  sight  normal  to  the  orbit. 

The  situation  with  the  dispersive  range  factor  W'R(y,  z)  of  (2.43)  is  somewhat  different.  The  maximum 
of  the  sine  is  attained  when  its  full  argument  is  zero,  i.e.,  when  Ry  =  cT(x°,  z,ui 0),  which  yields: 


Ry  —  R; 


=  cT(x°,  z,  w0)  -  Rz  =  Rz 


1  47re^ 

Wa  2  TOp 


■at. 


(2.44) 


This  is  obviously  not  equivalent  to  Ry  —  Rz  =  0,  which  implies  that  on  the  image  the  target  is  shifted  in 
range  from  its  absolute  position  by  the  value  of  the  right  hand  side  of  (2.44).  The  range  resolution  per  se, 
however,  is  still  given  by  the  width  of  the  main  lobe  of  the  sine  as  a  function  of  the  argument  Ry  —  Rz: 


A  R  = 


7 rc 
a'r 


27 rc  /  2St\ 


(2.45) 


The  factor  of  2  in  front  of  St  in  equation  (2.45)  appears  because  the  pulse  dilation  is  evaluated  for  the 
round-trip  between  the  antenna  and  the  target.  The  relative  deterioration  of  the  range  resolution  compared 
to  the  non-dispersive  case  is  therefore  ~  2 St/t,  see  (2.45).  As  shown  in  [Tsy09a],  for  the  typical  values  of 
all  the  relevant  parameters  in  the  deterministic  ionosphere  this  quantity  is  very  small,  about  0.02%. 

If  the  quadratic  term  is  retained  in  the  exponent  under  the  integral  (2.42)  after  the  change  of  variables, 
then  the  integration  yields  a  complex- valued  function  WR(y,z)  rather  than  the  actual  sine.  As  shown 
in  [Tsy09a],  the  maxima  and  minima  of  its  absolute  value  \  WR(y,  z)  |  are  located  precisely  where  the  maxima 
and  zeros  of  the  sine  are.  In  this  sense,  the  range  resolution  remains  unchanged.  However,  instead  of  the 
actual  zeros  the  new  \WR(y,  z)\  now  has  minima  with  some  nonzero  values.  In  particular,  for  the  first 
minimum  we  have 


\WR(y,z)\ 


ia'(Ry/c-T(x0,z,ui0))  =  ^- 


=  t~^2T  0  +  ^{SaT2))  , 


(2.46) 


The  right-hand  side  of  equality  (2.46)  is  about  0.02  •  r  for  the  typical  parameters  involved.  Therefore,  the 
dispersion  of  radio  waves  in  the  ionosphere  causes  a  deterioration  of  the  image  sharpness.  As  the  central 
maximum  of  \WR(y,  z)|  basically  remains  equal  to  r,  the  extent  of  deterioration  appears  2%,  but  this  value 
can,  of  course,  vary  both  ways  depending  on  the  parameters  of  the  radar  and  the  conditions  of  the  ionosphere. 

Let  us  also  recall  that  the  quantity  A f  that  determines  the  dilation  of  the  pulse  and  its  new  rate,  see 
formulae  (2.39)  and  (2.30),  as  well  as  the  shift  of  the  target  position,  see  formula  (2.44),  does  not  only 
contain  a  deterministic  component,  but  also  a  stochastic  one.  The  latter  is  the  average  value  of  y  along  the 
ray  trajectory,  see  (2.37).  Hence,  based  on  the  idea  of  ergodicity,  and  recalling  that  (y)  =  0,  we  can  expect 
that  the  effect  of  the  aforementioned  stochastic  component  on  the  radar  performance  in  range  will  not  be 
large. 

To  quantify  this  argument  further,  we  denote  the  random  contribution  to  TgI(x° ,  z,  w)  by 


R-. 


‘fio  =  J  y{s)ds , 


Rz  =  \x°-z\,  (ip0)  =  0, 


(2.47) 


where  the  integration  is  conducted  along  the  unperturbed  trajectory,  and  also  introduce  the  mean  dilation 
according  to  (2.39): 


(St) 


2 Rz  Ane2  B 
c  meu!Q  ojq 


<A0 


2 Rz  Aire2  B  N^H) 


c  meu>Q  wq  H 
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Standard  deviations  of  both  ipo  and  St  were  computed  in  [Tsy09a]  using  the  spectrum  of  turbulent  fluctua¬ 
tions  (2.11).  It  was  also  shown  that  > J (ip^)  <C  and  \J {St2)  <C  (St),  or  more  precisely,  for  the  case  of  a 
homogeneous  ionosphere, 


V^o) 

n^h) 


(St) 


1, 


where  tq  is  the  correlation  length,  and  M  is  given  by  (2.9).  In  the  case  of  a  stratified  ionosphere,  the  variance 
of  fluctuations  becomes  a  function  of  the  altitude  h:  (y2)  =  (y2(h)),  and  one  can  show  [Tsy09a]  that 


where  N2'^H^ 


H 

J  (Ne(h))2dh. 
0 


(2.48) 


According  to  [Arm05],  the  value  of  N2,('H'1  for  high  altitudes  H  is  between  9.3 •  1018cm-5  and  9.9  •  102Ocm-5, 
with  the  average  about  1.5  •  1019c?n-5.  Then,  say,  for  H  =  500 km  we  have  ■  H )1/2  «  2.7  •  1013cto-2, 

which  is  close  to  Ne  ~  (Ne)  ■  H  ss  5  •  1013c?n-2  if  one  takes  (Ne)  =  106cm-3.  We  therefore  conclude  that 
the  contribution  of  turbulent  fluctuations  to  travel  delays,  as  well  as  to  the  change  of  the  pulse  shape  during 
its  propagation  in  the  ionosphere,  is  much  smaller  than  that  of  the  baseline  dispersion.  Therefore,  the  role 
of  randomness  in  the  deterioration  of  SAR  performance  in  range  appears  minute,  and  can  be  disregarded  for 
the  typical  values  of  the  parameters  involved. 

2.4.2.  Azimuthal  resolution.  The  azimuthal  factor  (2.41)  of  the  generalized  ambiguity  function  has 
also  been  computed  in  [Tsy09a].  Similarly  to  (2.47),  let  us  denote 


n,z 

<Pn  =  J  y(s)ds,  Rz 


=  \X  —  z 


We  also  linearize  the  travel  distances: 


Rz  =  \z  -  xn\  =  \JRq  +  (zi  -  a:")2 
Ry  =  \y-  xn\  =  y/R%  +  (yi  -  x^)2 


(zD2 

R2 


iwn 

2  R2  )’ 


(y i  -  a:")2 

R2 


Ro 


1  {yi-xi)2\ 

2  Rl  J  ’ 


(2.49) 


(2.50) 


because  |a;"|  <C  Ro  and  \y\  —  x™\  <C  Ro,  where  Rq  is  the  distance  along  the  normal  (recall,  the  subscript  “1” 
denotes  the  coordinate  along  the  orbit,  and  Z\  =  0,  as  in  Section  2.3.1).  Then,  up  to  a  factor  of  magnitude 
one,  formula  (2.41)  along  with  (2.37b)  yields: 


2iWQ 

WA(y,z)  oc  J2  e 

n——N/2 


VI*?  |  l 
Ro  ^ 2 


(®?r 


*(H> 


2  Ro 


Wr 


(2.51) 


According  to  the  central  limit  theorem,  for  long  integration  distances  Rz  each  quantity  <pn  of  (2.49)  becomes 
a  Gaussian  random  variable  with  zero  mean,  see  [RKT89b,  Ch.  I],  Hence,  the  mean  of  each  term  of  (2.51)  is 


VI  *i  ,  1  (®i*)2  4-ire2  Ne. 11 1  i  1 
J?0  ^2  2 B0  ra^  '  2 
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(2.52) 
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and  the  variance  is  given  by  (see  [Tsy09a]): 


a 


2 

n 


Wn 


(—1) 


2 


{&)■ 


(2.53) 


The  (dimensionless)  quantity  on  the  right-hand  side  of  (2.53),  which  is  obtained  by  Taylor’s  formula,  is  the 
variance  of  the  phase;  it  is  small  for  the  typical  values  of  the  parameters  involved.  Thus,  the  decrease  of  the 
amplitude  on  the  right-hand  side  of  (2.52)  is  also  small;  in  the  literature  it  is  known  as  extinction  [RKT89b, 
Clr.  I], 

To  obtain  the  mean  and  variance  of  the  sum  (2.51),  we  will  assume  that  all  individual  variances  (2.53) 
are  the  same,  cr2  =  ct2,  n  =  — TV/2, . . .  ,N/2.  To  that  effect  we  note  that  while  the  integration  distance  Rz 
in  formula  (2.49)  depends  on  n,  its  variation  compared  to  R0  is  small,  see  (2.50). 

For  the  mean  value  of  (2.51)  we  can  therefore  write  using  (2.52): 


W(»,*)> 


oc 


n=-N/  2 


(  »1*1  |  1  Gl)2  4,e2 

V  R0  ~1"2  2Ro  H  J 


(2.54) 


where  the  deterministic  sum  on  the  right-hand  side  of  (2.54)  controls  the  azimuthal  resolution  in  the  case 
of  non-fluctuating  ionosphere.  There  are  two  terms  in  the  exponent  on  the  right-hand  side  of  (2.54).  The 
first  one,  linear  with  respect  to  x/,  is  the  only  thing  that  one  would  have  obtained  if  the  propagation  was 
unobstructed.  The  second  term,  which  is  quadratic  with  respect  to  a:",  is  due  to  the  phase  mismatch  between 
the  actual  received  signal  subject  to  ionospheric  delays  and  the  non-corrected  matched  filter. 

With  no  quadratic  term  we  have 
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n——N/2 


2iu>Q 

c  Rq 


( TV  +  l)sinc 


upyiAx^N+1) 

Roc 


(2.55) 


where  Aaq  is  the  distance  that  the  satellite  travels  along  the  orbit  between  the  successive  emission  of  pulses, 
x"  =  n  ■  Axi,  and  we  have  assumed  that  uloVI/^Xl  -C  1,  which  holds  for  the  typical  values  of  all  parameters. 
The  semi-width  of  the  main  lobe  of  the  sine  corresponds  to  the  value  of  its  argument  equal  to  7 r: 

2nyiAxi  2Af?o 
Rq  A  Ax  1 L 


where  we  have  replaced  TV  +  1  by  TV  = 


2XRq 

Ax±L 


2XRq 
Ax  1 L  ’ 


4?n  for  the  satellite  travel 


-,  because  TV  1  (Axi 

speed  of  8km/ sec).  Consequently,  the  azimuthal  resolution  in  the  undistorted  case  is  2yi  =  L/2. 

With  the  quadratic  term  taken  into  account,  the  function  ( W'A(y ,  z))  of  (2.54)  becomes  complex- valued. 
The  central  maximum  of  its  modulus  is  attained  when  the  argument  is  equal  to  zero,  and  the  first  minimum 
appears  precisely  where  the  first  zero  of  the  sine  (2.55)  is.  In  this  sense,  the  azimuthal  resolution  per  se  re¬ 
mains  unaffected.  However,  similarly  to  the  case  of  range  resolution,  the  sharpness  of  the  image  deteriorates, 
and  this  deterioration  for  the  azimuthal  resolution  appears  stronger.  Specifically,  one  can  show  [Tsy09a]  that 


\(w'A(y,z))\ 


u>0y1Ax1(N  +  l) _ 

ROc 


wo(Axi)2  47re2  TVi^  (TV  +  l)2 

2R0c  TOeWp  H  2ir2 


(2.56) 


which  is  approximately  0.09  •  (TV  +  1)  for  the  typical  values  of  the  parameters  involved.  Hence,  there  may  be 
about  9%  degradation  of  image  sharpness  in  the  azimuthal  direction  due  to  the  dispersion  of  radio  waves  in 
the  ionosphere. 

The  effect  of  randomness  on  the  azimuthal  resolution  is  also  stronger  than  that  on  the  range  resolution. 
Recall,  range  resolution  is  basically  determined  by  a  single  look,  and  because  of  the  ergodicity  the  contribution 
of  turbulent  fluctuations  into  the  radar  reading  on  this  single  look  effectively  vanishes  (i.e.,  averages  out).  On 
the  contrary,  azimuthal  resolution  is  determined  based  on  multiple  looks,  n  =  —TV/2, . . . ,  TV/2,  and  whereas 
the  contribution  of  randomness  into  each  individual  reading  is  still  small,  it  can  get  amplified  when  those 
individual  readings  are  summed  up.  In  [Tsy09a,  ST10],  we  call  this  effect  “un-averaging.” 


16 


S.  V.  Tsynkov 


To  estimate  the  variance  of  the  sum  (2.51)  one  needs  to  see  to  what  extent  the  individual  random  variables 
that  compose  this  sum  are  dependent  or  independent.  Each  of  these  random  variables  has  a  logarithmically 
normal  probability  distribution,  and  similarly  to  the  standard  Gaussian  variables,  if  they  are  uncorrelated, 
they  are  independent  [Tsy09a,  App.  F].  Therefore,  the  covariance  (correlation  function)  of  the  received  held 
along  the  synthetic  antenna  can  provide  a  measure  of  dependence  or  independence  for  the  random  variables 
that  compose  the  sum  (2.51).  In  [RKT89b,  Ch.  I]  it  is  shown  that  if  the  fluctuations  of  the  phase  are 
small,  then  the  correlation  function  for  the  held  approximately  coincides  with  the  correlation  function  for 
the  electron  number  density  (2.8).  Consequently,  the  correlation  length  ro,  see  formula  formula  (2.11),  will 
also  provide  a  characteristic  scale  of  how  rapidly  the  received  held  decorrelates  along  the  synthetic  antenna. 
Roughly  speaking,  for  the  locations  that  are  further  apart  than  r0  the  received  pulses  will  be  uncorrelated 
and  hence  independent,  whereas  for  the  locations  that  are  closer  than  ro  they  will  not  be  independent. 

Based  on  these  considerations,  the  following  expression  was  obtained  in  [Tsy09a]  for  the  variance  of 
(2.51): 

4  =  -N2ca2  =  ^(N+  1 )  V,  (2.57) 

ro  D 

where  D  =  Rq^-  is  the  length  of  the  synthetic  array.  It  has  also  been  shown  in  [Tsy09a]  that  the  correspond¬ 
ing  relative  change  of  the  azimuthal  resolution  is  |  ■\f%-  This  quantity  is  about  0.3%  for  the  typical  values 

of  the  parameters,  and  it  is  also  proportional  tow0  ,  i.e. ,  increases  as  the  carrier  frequency  decreases. 

A  0.3%  degradation  of  the  azimuthal  resolution  is  not  negligible  but  is  still  rather  small.  However,  when 
estimating  cr  via  (<p2),  see  formulae  (2.48)  and  (2.53),  we  have  taken  the  average  values  for  M  and  N^H- 
Taking  the  corresponding  values  from  the  upper  part  of  their  respective  ranges  may  result  in  a  much  larger 
value  of  (ip2)  and  as  such,  in  a  noticeable  deterioration  of  azimuthal  resolution. 

2.5.  Dual  carrier  probing.  To  remove  the  mismatch  that  causes  distortions  of  the  image,  one  needs 
to  correct  the  filter,  i.e.,  replace  the  unobstructed  travel  times  | y  —  xn\/c  by  Tgr(xn,  y,  w0)  in  (2.40)  and 
by  Tph(a:",  y,u>o)  hr  (2-41),  and  also  adjust  the  chirp  duration  and  chirp  rate  in  the  filter  factor  under  the 
integral  (2.40).  However,  unlike  in  the  received  field,  which  is  a  physical  observable,  the  correction  in  the 
filter  must  be  done  theoretically.  Therefore,  one  needs  to  know  the  quantity  A f,  see  formula  (2.38),  that 
characterizes  the  ionosphere.  The  availability  of  J\f  will  allow  one  to  calculate  the  dilation  (2.39)  and  the 
new  chirp  rate  (2.30)  for  any  reference  location  y. 

The  quantity  N  depends  on  the  deterministic  part  of  the  electron  number  density,  see  formula  (C.8), 
and  also  has  a  random  component  obtained  by  integrating  y(x)  along  the  ray  trajectory.  In  the  literature, 
many  estimates  are  available  for  the  electron  number  density  in  the  ionosphere,  see,  e.g.,  [Gin64]  or  [Bud85]. 
Those  estimates,  however,  can  only  provide  a  typical  range  of  values,  especially  as  the  parameters  of  the 
ionosphere  are  known  to  vary  in  space  and  in  time.  This  does  not  allow  one  to  correct  the  filter  with  a 
sufficient  degree  of  reliability.  What  is  rather  needed  is  an  accurate  value  of  AT  exactly  at  the  time  and  place 
the  image  is  taken. 

Hereafter,  we  exploit  the  idea  of  dual  carrier  probing  for  deriving  Af .  Let  us  assume  that  there  is  an 
object  or  feature  in  the  scene  that  can  be  clearly  identified  on  the  image.  This  object  does  not  have  to  be 
artificial.  It  does  not  have  to  dominate  the  scene,  say,  by  having  the  highest  reflectivity.  Its  location  does 
not  have  to  be  known  ahead  of  time.  It  merely  has  to  be  something  that  can  be  fairly  easily  picked  out  and 
matched  on  different  images  that  represent  the  same  terrain.  For  example,  it  can  be  some  landmark,  such 
as  a  hilltop,  a  building,  a  road  intersection,  etc. 

Let  wo  and  wi  be  two  distinct  carrier  frequencies,  wo  ^  u>i,  and  let  and  R^  be  the  corresponding 
ranges  of  the  aforementioned  object  measured  by  the  radar,  whereas  its  true  range  is  Rz  (unknown  yet). 
Then,  we  can  consider  the  corresponding  two  equations  (2.44)  as  a  system: 

=cTgr(£c°,z,w0), 

4  ’  (2-58) 

Ry]  =cTgI{x°,z,cj1), 

where  the  group  travel  time  Tgr  is  given  by  formula  (2.37a).  With  the  data  R^  and  Ay  1  available,  system 
(2.58)  can  be  solved  with  respect  to  the  two  unknown  quantities:  the  true  range  Rz  and  the  integral  quantity 
AT  that  characterizes  the  plasma. 


AFOSR  Grant  No.  FA9550-07-1-0170.  Final  Performance  Report 


17 


More  precisely,  the  true  range  Rz  enters  into  the  expression  for  Tgr  as  a  factor  in  front  of  the  formula, 
and  also  as  the  integration  limit  for  the  random  component  of  the  electron  number  density,  see  (2.37a). 
As,  however,  we  have  seen  in  Section  2.4,  the  contribution  of  randomness  into  the  single-look  image  is  very 
small  for  large  Rz  because  of  the  ergodicity.  Randomness  rather  manifests  itself  when  the  information  from 
multiple  looks  is  combined  for  the  analysis  of  azimuthal  resolution.  Therefore,  since  system  (2.58)  contains 
only  single-look  information,  we  can  disregard  the  stochastic  contribution  into  Tgr  and  instead  of  (2.37a) 
substitute  the  expression 


TgI (x°,  z,  u) 


Rz  I-  1  47re2  jVjg) ' 

c  [  uj2  2 me  H 


(2.59) 


on  the  right  hand  side  of  equations  (2.58).  We  emphasize  that  the  simplification  offered  by  formula  (2.59) 
compared  to  formula  (2.37a)  is  due  to  the  ergodicity  of  y  that  can  be  employed  in  the  context  of  single-look 

_  /  rr\ 

imaging.  This  simplification  is  important  as  it  allows  one  to  determine  the  plasma  quantity  Ne  ;  defined  by 
formula  (C.8),  as  well  as  the  true  range  Rz. 

Let  us  introduce  the  following  notation  for  brevity: 


wpe  = 


47re2  Ne H'> 


H 


(2.60) 


Using  this  notation,  system  (2.58)  with  the  expression  for  Tgr  given  by  (2.59)  is  recast  as 


(i+S)- 


(2.61) 


The  range  Rz  can  be  eliminated  from  system  (2.61)  by  dividing  the  equations  by  one  another.  Then,  after 
simple  transformations,  we  have: 


2^0V(<)  -  41}) 


(2.62) 


We  note  that  the  denominator  in  formula  (2.62)  does  not  turn  into  zero,  because  if  it  was  equal  to  zero,  then 

_  /  rr\ 

system  (2.61)  would  immediately  yield  u) o  =  The  actual  electron  content  Ne  ;  can  be  determined  using 
(2.60)  once  w2e  has  been  evaluated  according  to  (2.62).  Moreover,  if  there  is  a  need  to  know  the  true  range 
Rz  of  the  chosen  reference  object,  it  can  be  easily  found  from  any  of  the  equations  (2.61). 

2.6.  Correcting  the  matched  filter.  Once  the  quantity  A'(:,/  has  been  found,  one  can  use  formu¬ 
lae  (2.37)  and  accurately  evaluate  the  deterministic  contributions  to  the  travel  times  Tph(xn,  y,  w0)  and 
TgT(xn,  y,u o)  for  any  location  y.  Then,  one  can  also  compute  the  pulse  dilation  for  the  reference  point  y 
using  formula  (2.39),  and  the  new  pulse  rate  according  to  (2.30).  Those  quantities  are  needed  to  modify  the 
definition  of  the  matched  filter  in  formulae  (2.40)  and  (2.41),  i.e.,  to  correct  the  matched  filter. 

2.6.1.  Range  resolution.  The  matched  filter  is  corrected  by  using  the  actual  travel  times  instead  of 
the  travel  times  in  vacuum,  and  by  taking  into  account  the  dilation  and  the  rate  change  of  the  pulse  upon 
its  arrival.  With  this  correction  in  place,  the  range  factor  of  the  generalized  ambiguity  function  becomes  [cf. 
formula  (2.42)]: 


W'R(y,z)  = 


oc 


rmin{T'(2/)/2+2fgr(2/),  r' (z)/2+2Tgr (z)} 

/max{  —  T'(y)/2+2fgr(y),  -r'(z)/2+2Tgr(z)} 
rmin{r'  (y)/2+2fgI(y),  t'  (z)/2+2Tgr  (z)} 


/max{  —  T,(y)/2+2Tei(y),  -r'(z)/2+2Tgr(z)} 


e— ia'(y)(t- 2TgI(y))2  eioc'(z)(t- 2Tgr(z))2^ 
ei(a'(z)-a'(j/))t2e4i(a'(j/)fgr(j/)-a'(z)Tgr(z))^^ 


(2.63) 
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where  we  have  introduced  the  abbreviated  notations  Tgr(z)  =  Tgl  (x° ,  z,uo)  and  TgI(y)  =  Tgr(x°,  y,  u>o). 
We  emphasize  that  for  the  filter  part  of  the  generalized  ambiguity  function  we  can  substitute  only  the 
deterministic  component  of  the  group  travel  time  (2.59),  because  this  is  what  we  derive  from  the  single-look 
information,  see  Section  2.5.  The  primed  quantities  t'  and  a'  in  formula  (2.63)  are  to  be  evaluated  for  the 
pulse  round-trip  between  the  antenna  and  the  target  [cf.  formulae  (2.29)  and  (2.30)]: 

t'  =  t  +  2 St  and  a'  =  a  —  St—t  =  —  f  1 - ]  ,  (2.64) 

r-  2 r  V  t  / 


where  Sr  is  given  by  formula  (2.39).  For  the  integral  (2.63),  we  consider  the  integration  limits  ^fr'(y)/2 
2 T(y),  and,  changing  the  integration  variable:  u  =  t  —  2 T(y),  obtain: 


WR{y,z)  oc 


/ 

/-r'(»)/ 2 
rT'(l 0/2 

/ 

'~r'(y)/2 


z-ioc'{y)u2  e*a'(z)(«+2fgr(y)-2Tgr(z))2 


,i(a'  (z)-a'  (y))u2  e4iot'  (z)(fgI(y)-Tgr(z))u 


(2.65) 


where  the  constant  factors  of  magnitude  one  in  front  of  the  integrals  are  dropped. 

Following  the  analysis  of  Section  2.4.1,  we  first  disregard  the  quadratic  term  ~  u 2  in  the  exponent  under 
the  integral  (2.65),  because  this  term  is  small.  Then,  integral  (2.65)  evaluates  to 

W'R{y ,  z)  oc  t  (y)smc[2a' (z)(Tgl(y)  -  TgI{z))r' (y)}.  (2.66) 


The  range  resolution  of  the  radar  is  defined  as  the  full  width  of  the  sine  function  (2.66)  interpreted  as  a 
function  of  Ry  —  Rz,  and  hence  can  be  obtained  by  setting  the  argument  of  the  sine  equal  to  27t: 


2o/(z)(fer(v)-Tg'(z))T'(y)=2n. 


With  the  help  of  (2.64),  the  previous  equality  transforms  into 


B 

T 


(Tgr{y)  -Tgr(z))r 


=  27T. 


(2.67) 


We  first  notice  that  the  correction  due  to  the  dilation  of  the  pulse  in  formula  (2.67)  becomes  quadratic  with 
respect  to  — .  For  the  typical  values  of  the  parameters  involved,  we  have  —  ~  10-4,  see  [Tsy09a,  ST10]. 
Consequently,  we  can  disregard  this  correction  completely  and  instead  of  (2.67)  write: 

B(fgI(y)-TgI(z))  =  27r.  (2.68) 

Then,  we  can  use  formulae  (2.37a),  (2.47),  (2.59),  and  (2.60),  and  express  the  difference  between  the  group 
travel  times  in  (2.68)  as  follows: 


TM  -  Tgr(Z) 


1  47re2  1 

2  meu> q  c 


<Po(z). 


(2.69) 


As  {(po(z))  =  0,  we  can  say  that  for  every  particular  realization  of  the  random  held  y  we  have  |<^o(-z))|  ~ 
\J (v2q).  Hence,  the  second  term  on  the  right-hand  side  of  equality  (2.69)  appears  much  smaller  than  the 
first  term,  because  the  exact  same  argument  applies  here  as  given  in  Section  2.4.1,  see  the  discussion  around 
equation  (2.48).  Thus,  disregarding  this  second  term  and  substituting  the  result  into  (2.68),  we  arrive  at  the 
following  expression  for  the  range  resolution: 


A R  =  Ry  -  Rz 


27 rc  /  _  1  ^pc\ 

~B~{  "2  ^4)' 


(2.70) 


This  is  an  improvement  over  (2.45)  and  even  a  slight  improvement  over  the  non-dispersive  resolution  A R  = 
2irc/B,  because  the  center  frequency  group  velocity  that  essentially  enters  into  the  numerator  of  formula 
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(2.70)  is  slower  than  the  speed  of  light  c.  However,  as  indicated  in  Section  2.4.1,  the  range  resolution  of  a 
SAR  sensor  does  not  suffer  much  from  the  ionospheric  dispersion  in  any  event.  Therefore,  what  is  even  more 
important  is  that  when  the  filter  is  corrected,  the  target  is  no  longer  shifted  in  range  from  its  true  position 
as  in  formula  (2.44).  That  is  because  the  maximum  of  the  sine  in  formula  (2.66)  is  attained  precisely  at 
Ry  =  Rz- 

What  is  also  very  important  is  that  the  degradation  of  image  sharpness  becomes  negligible  once  the 
filter  has  been  corrected.  To  analyze  this  effect,  we  bring  back  the  quadratic  term  in  the  exponent  under  the 
integral  (2.65).  First,  we  recast  this  integral  in  the  form 


where 


rT'(y)/2 

W'R(y,z)  oc  /  eilu  eKudu, 

J-T'{y)/2 


7  =  a'{z)  -  a'{y) 


^(5r(y)  -  6t(z)) 


B  2 (Ry  -  Rz)  47re2  B 

T 2  C  ?7leWQ  W 0 


(2.71) 


(2.72) 


C  —  4a  ( z)(TgI(y )  Tgr(z)), 


(2.73) 


and  the  quantity  A f  is  defined  in  formula  (2.38).  The  quadratic  term  in  the  exponent  under  the  integral 
(2.71)  is  small:  |7W2|  <C  1  for  u  €  [— r'(y)/2, r'(j/)/2],  and  consequently  we  have: 


Wjiiy.z) 


'(!/)/ 2 


l-r'(y)/ 2 


(1  +  i'yu2)el^udu 


=  r'  sine  [2  c/  ( z ) (Tgr (y)  -  Tgr(z))r']  +  *7 


.  4Cr'cos^ 


(-8  +  CV2)sin^ 


(2.74) 


where  t'  =  r'(y).  The  first  term  on  the  right-hand  side  of  (2.74)  obviously  coincides  with  (2.66),  and  the 
second  term  represents  a  correction.  According  to  (2.72)  and  (2.73),  (2.69),  both  7  and  /  vanish  as  Ry 
approaches  Rz.  However,  the  fraction  on  the  right-hand  side  of  (2.74)  remains  bounded  as  (  — >  0  (and 
7  — >  0).  It  is,  in  fact,  easy  to  see  that 

W'R{y,  z)  «  r'sinc [2c/ ( z ) (Tgr (y)  -  Tgr(z))r']  +  (2.75) 

As  the  first  term  on  the  right-hand  side  of  (2.75)  is  proportional  to  r',  the  relative  magnitude  of  the  correction 
is  ~  yr'  (this  is  a  dimensionless  quantity).  There  is  no  correction  at  the  central  peak  of  the  sine,  because 
7  =  0  when  Ry  =  Rz.  To  quantify  the  extent  to  what  the  image  sharpness  is  affected,  we  estimate  the 
correction  at  the  first  zero  of  the  sine.  Using  formulae  (2.72)  and  (2.64),  we  obtain: 


7T,2  =  B2 

c  mew 0  W0  V  T  ) 

where  the  value  of  Ry  —  Rz  shall  be  taken  according  to  formula  (2.70).  For  the  typical  values  of  the  parameters 
involved,  we  have 


47re2  cu2 

O  ~  9 

meuj  0  wg 


10  4  and  — 
Wo 


10 


-2 


so  that  altogether  we  can  write: 


3 

16 


7  T 


37T  w|e  B_  ( -.  _  1  f,  4St 
4  w2  w0  ^  2  w2  J  (  +  r 


37 r 

T 


10"6. 


We  therefore  conclude  that  the  degradation  of  image  contrast  in  range  decreases  to  approximately  0.0002% 
as  opposed  to  2%  in  the  non-corrected  case,  see  formula  (2.46). 
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2.6.2.  Azimuthal  resolution.  Once  the  matched  filter  has  been  corrected,  expression  (2.41)  for  the 
azimuthal  component  of  the  generalized  ambiguity  function  becomes 


N/2 


W'A(y,  z)  =  J2 


o2iu0(Tph(xn  ,y,uo)  —  Tph(xn,z,u0)) 


(2.76) 


n——N/2 

where  in  the  filter  part  we  have  substituted  the  deterministic  portion  of  the  phase  travel  time: 


Tph(x°,z,uj)  =  — 
c 


1  - 


1  47re2  NgH^ 
lo2  2me  H 


because  this  is  what  we  derive  from  the  sing-look  image  by  means  of  dual  carrier  probing,  see  Section  2.5. 
Then,  using  formulae  (2.37b),  (2.49),  and  (2.50),  we  can  transform  expression  (2.76)  into  [cf.  formula  (2.51]: 


W'A(y,z)  oc  E  e 

n=-N/2 


N/2  2iui Q  |  _  Bl^j* 

HO 


1-1- 


+  9  - 


~^ipr 


(2.77) 


where  the  constant  factor  of  unit  magnitude  in  front  of  the  sum  was  dropped.  Next,  following  the  same 
derivation  as  in  Section  2.4.2  (see  also  [Tsy09a]),  we  compute  the  mean  and  variance  of  each  term  in  the 
sum  (2.77): 


i-l- 


AH) 


+  9  - 


(2.78) 


(vl) 


<7 


2 

n 


U  n 


(1^) 


(vl)- 


(2.79) 


Also  as  in  Section  2.4.2,  the  individual  variances  (2.79)  can  be  assumed  approximately  equal  to  one  another, 
<r2  s=s  a2,  n  =  —N/2, . . . ,  N/2,  because  the  integration  distance  for  ipn  depends  on  n  weakly.  Hence,  we  can 
use  (2.78)  and  derive  the  mean  value  of  the  sum  (2.77): 


(WA(y,z)) 


(X 


N/2 

e ~°2/2  E 


n=  —  N/2 


-AH) 


)]. 


(2.80) 


The  sum  (2.80)  can  be  actually  computed  using  the  same  approach  as  employed  in  Section  2.4.2  and  in 
[Tsy09a].  The  result  is  another  sine  function,  and  it  leads  to  the  following  estimate  of  the  azimuthal  resolution 
in  the  case  of  a  non-fluctuating  ionosphere: 


L  (  1  47re2  NlH)\ 

2  y  2  meWQ  H  J 


This  value  is  only  marginally  worse  than  that  obtained  in  the  non-dispersive  case:  2yi  =  L/2.  Even  more 
important,  correction  of  the  filter  restores  the  sharpness  of  the  image  for  the  non-fluctuating  ionosphere,  so 
that  there  is  no  deterioration  like  in  the  non-corrected  case,  see  formula  (2.56).  The  reason  is  that  unlike  in 
formula  (2.51),  there  is  no  quadratic  term  with  respect  to  a;"  in  the  exponent  in  the  sum  (2.77). 

As,  however,  the  proposed  correction  is  based  on  the  single  look  data  that  essentially  do  not  account  for 
randomness  (due  to  the  ergodicity) ,  it  does  not  offer  any  remedy  for  the  stochastic  part  of  image  distortions 
in  the  azimuthal  direction.  In  other  words,  we  can  expect  that  even  with  the  corrected  filter  in  place,  the 
azimuthal  factor  WA(y,  z)  of  the  generalized  ambiguity  function  will  still  have  the  variance  of  order  (2.57), 
which  will  result  in  the  same  extent  of  deterioration  in  the  azimuthal  resolution  as  estimated  in  [Tsy09a]. 
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2.7.  Discussion.  SAR  images  obtained  from  satellites  are  prone  to  deterioration  due  to  the  temporal 
dispersion  of  radio  waves  in  the  Earth’s  ionosphere.  The  deterioration  is  stronger  for  lower  carrier  frequencies 
and  weaker  for  higher  carrier  frequencies.  We  have  analyzed  this  phenomenon  in  the  case  of  a  scalar  field 
(polarization  of  radar  pulses  is  not  taken  into  account)  propagating  in  an  inhomogeneous  cold  plasma.  Our 
analysis  shows  that  image  deterioration  is  due  to  the  mismatch  between  the  actual  signals  scattered  off  the 
Earths’  surface  and  received  by  the  radar  antenna  and  the  matched  filter  (an  image  processing  tool)  designed 
as  if  the  propagation  between  the  antenna  and  the  ground  was  unobstructed. 

To  correct  the  filter,  one  needs  to  know  a  key  characteristic  of  the  ionosphere  precisely  at  the  time  and 
place  the  image  is  taken.  This  key  characteristic  is  the  integral  content  of  electrons  per  unit  area  in  the  layer 
of  ionospheric  plasma  between  the  satellite  and  the  ground.  We  have  proposed  to  derive  this  quantity  by 
probing  the  terrain,  and  hence  the  ionosphere,  on  two  distinct  carrier  frequencies.  We  have  also  shown  that 
the  resulting  correction  of  the  filter  completely  eliminates  all  the  distortions  of  the  image  that  are  due  to  the 
deterministic  part  of  the  electron  content  in  the  ionosphere. 

The  correction,  however,  is  not  efficient  for  removing  the  distortions  due  to  the  random  part  of  the 
charged  particle  content,  i.e. ,  turbulent  fluctuations  of  the  electron  number  density.  The  reason  is  that  the 
correction  is  based  on  the  data  from  a  single-look  image,  for  which  the  effect  of  randomness  is  “downplayed” 
by  ergodicity.  Mitigating  the  stochastic  part  of  image  distortions  will  be  a  subject  for  the  future  study. 

In  addition,  we  should  note  that  our  approach  to  correcting  the  filter  relies  on  an  important  consideration 
that  the  mean  parameters  of  the  ionosphere  do  not  vary  in  the  horizontal  direction.  This  assumption  will 
hold  only  if  the  size  of  the  imaged  scene  is  relatively  small.  For  larger  scenes,  the  procedure  of  dual  carrier 
probing  needs  to  be  carried  out  repeatedly,  and  the  filter  needs  to  be  adjusted  accordingly,  as  the  satellite 
travels  along  the  orbit  and  the  antenna’s  footprint  sweeps  the  imaged  stripmap. 

The  issues  to  be  studied  in  the  future  may  include: 

•  An  investigation  of  whether  the  stochastic  component  of  image  distortions  can  be  tackled  by  any 
of  the  “black  box”  image  sharpening  and  deblurring  techniques  developed  previously  in  the  general 
framework,  see,  e.g.,  [GW08],  with  no  direct  relation  to  spaceborne  SAR.  The  considerations  based 
on  the  central  limit  theorem  that  are  outlined  in  Section  2.4.2  may  be  useful  in  this  perspective. 

•  Introduction  of  a  vector  propagation  model,  analysis  of  SAR  imaging  with  polarization,  accounting 
for  anisotropy  of  the  medium  due  to  the  magnetic  field  of  the  Earth,  and  for  the  Faraday  rotation. 

•  Analysis  of  the  effect  of  turbulent  fluctuations  in  the  ionosphere  on  the  image  using  true  Kolmogorov- 
type  spectra,  as  opposed  to  the  modified  spectrum  (2.11)  that  accounts  only  for  the  short  range 
phenomena.  Accounting  for  anisotropic  turbulence  due  to  the  magnetic  field  of  the  Earth. 

•  Approaches  for  reducing  the  noise  that  is  always  present,  regardless  of  whether  the  imaging  is  done 
through  the  ionosphere  or  not,  for  example,  instrument  noise  or  rough  terrain  noise.  As  a  minimum 
requirement,  this  noise  may  not  be  amplified  by  any  of  the  techniques  used  for  mitigating  the 
ionospheric  distortions. 

•  Analysis  of  the  dispersion  on  the  target,  and  distinguishing  between  the  target  dispersion  and  the 
ionospheric  dispersion.  Possible  applications  to  material  identification  SAR  (MISAR). 

•  Optimization  of  interrogating  waveforms  for  the  ionosphere,  beyond  the  standard  linear  upchirps. 

3.  On  the  use  of  start-stop  approximation  for  spaceborne  SAR  imaging. 

3.1.  Background.  A  synthetic  aperture  radar  (SAR)  builds  an  image  of  a  target  (could  be  an  area 
of  the  Earth’s  surface)  by  successively  illuminating  it  with  a  series  of  electromagnetic  pulses  transmitted 
from  different  locations,  and  subsequently  processing  the  scattered  waves  with  the  help  of  a  matched  filter. 
Typically,  the  antenna  of  such  a  radar  is  mounted  on  an  aircraft  (airborne  radar)  or  a  satellite  (spaceborne 
radar),  and  the  interrogating  pulses  are  emitted  and  received  as  the  craft  travels  along  its  flight  track  or 
orbit,  respectively.  A  standard  approach  employed  when  processing  the  SAR  data  is  based  on  the  start-stop 
approximation,  when  the  antenna  is  assumed  at  standstill  while  it  sends  a  pulse  and  receives  the  scattered 
response,  after  which  it  moves  on  to  the  position  where  the  next  pulse  is  emitted,  see,  e.g.,  [CM91,CGM95, 
FL99]  and  also  [CheOl].  This  approximation  is  intuitively  well  justified  for  an  airborne  radar.  Indeed,  a 
typical  aircraft  speed  is  between  200?n/sec  and  300?n/sec.  Then,  even  when  the  look  angle  is  large  (close  to 
7r/2)  so  that  the  radar  scans  the  surface  far  to  the  side  rather  than  directly  underneath  the  flight  track,  the 
round-trip  travel  time  for  the  pulse  will  still  be  under  10_3sec.  This  yields  the  displacement  of  the  antenna 
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between  20cm  and  30c?n,  see  Section  3.3,  which  is  on  the  order  of  one  to  a  few  wavelengths  for  the  microwave 
frequency  band  that  the  radar  typically  operates  in,  and  is  also  far  shorter  than  the  typical  antenna  size. 
Besides,  the  magnitude  of  the  linear  Doppler  frequency  shift  in  this  case  is  rather  small,  under  10  ,  see 

Section  3.4. 

On  the  other  hand,  using  the  start-stop  approximation  for  a  spaceborne  radar  may  raise  questions.  With 
the  round-trip  travel  times  ~  l/200sec  and  the  speed  of  8 km/ sec,  the  displacement  of  the  antenna  over  the 
period  of  time  between  the  emission  and  arrival  of  a  pulse  may  reach  about  40?n,  which  is  several  times 
longer  than  the  antenna  size  and  much  larger  than  the  wavelength.  It  is  not  clear  ahead  of  time  whether 
neglecting  this  displacement  may  affect  the  quality  of  the  image.  Moreover,  the  Doppler  frequency  shift, 
which  is  ~  v/c,  where  v  is  the  speed  of  the  satellite  and  c  is  the  speed  of  light,  see  Section  3.4,  will  also  be 
about  one  and  a  half  orders  of  magnitude  larger  than  in  the  case  of  an  airborne  radar.  In  the  earlier  SAR 
studies,  the  Doppler  shift  was  used  to  explain  the  mechanism  of  along-the-track  resolution,  see,  e.g.,  [FL99]. 
However,  when  building  a  generalized  ambiguity  function  the  Doppler  effect  is  often  neglected,  and  it  is  not 
apparent  whether  it  may  cause  a  deterioration  of  the  image. 

In  this  part  of  the  report,  we  analyze  the  generalized  ambiguity  function  of  a  SAR  system  while  taking 
the  two  foregoing  phenomena  into  account.  To  enable  the  analysis,  we  need  to  re-derive  the  non-dispersive 
generalized  ambiguity  function  of  Section  2.3.1  for  the  case  when  the  sending  and  receiving  locations  can 
be  different,  see  Section  3.2.  Our  analysis  subsequently  shows  that  the  displacement  of  the  antenna,  even 
though  quite  noticeable,  has  a  very  small  detrimental  effect  0(y/c )  on  the  resolution  of  the  image  and  can 
therefore  be  disregarded  for  all  practical  purposes,  see  Section  3.3.  As  for  the  Doppler  effect,  disregarding 
it  in  the  generalized  ambiguity  function  may  cause  an  overall  translation  of  the  entire  scene  on  the  image, 
as  well  as  a  slight  deterioration  of  the  sharpness  (contrast),  see  Section  3.4.  To  avoid  this,  one  must  include 
the  frequency  shift  into  the  definition  of  a  matched  filter. 

3.2.  The  generalized  ambiguity  function.  For  the  analysis  in  this  section,  we  will  employ  the 
following  assumptions.  The  electromagnetic  field  propagating  between  the  antenna  and  the  target  will  be 
assumed  scalar  and  will  be  governed  by  the  standard  d’Alembert  equation  with  the  speed  c  (the  speed  of 
light).  In  other  words,  we  will  disregard  the  polarization.  The  image  will  be  assumed  two-dimensional, 
i.e.,  the  radar  will  measure  coordinates  of  the  targets  along  the  Earth’s  surface,  but  will  not  measure  their 
elevation.  All  targets  will  be  assumed  dispersionless,  and  scattering  off  the  targets  will  be  linearized  and 
analyzed  using  the  first  Born  approximation,  see  [BW99,  Section  13.1.4]. 

The  radar  will  be  assumed  to  operate  in  the  standard  mode  known  as  stripmap,  see  [CGM95].  In  this 
case,  the  antenna  points  in  the  direction  which  is  fixed  relative  to  that  of  the  satellite  motion.  Hence,  the 
footprint  of  the  beam  emitted  by  the  antenna  sweeps  a  strip  on  the  Earth’s  surface  parallel  to  flight  track,  i.e., 
to  the  orbit.  In  doing  so,  the  antenna  does  not  necessarily  have  to  point  precisely  sideways,  i.e.,  normally  to 
the  direction  of  motion.  Instead,  it  can  point  either  forward  or  backward  at  a  fixed  angle,  which  corresponds 
to  the  so-called  squinted  stripmap  mode. 

Finally,  we  will  need  to  use  the  Lorentz  transform  [Gar07,  Chapter  3,  Section  3]  that  preserves  the  form 
of  the  d’Alembert  equation  in  the  case  of  moving  transmitters  and/or  receivers,  see  Section  3.4.  To  be  able 
to  do  so,  we  will  have  to  assume  that  the  motion  of  the  antenna  is  straightforward  and  uniform.  Technically 
speaking,  this  is  not  true  for  a  satellite  orbiting  the  Earth.  However,  for  a  relatively  short  stretch  of  the 
orbit,  along  which  the  measurements  are  taken  of  a  given  small  target  (i.e.,  for  which  a  given  target  stays  in 
the  beam,  see  Section  3.3),  this  is  a  sufficiently  accurate  approximation. 

Hereafter,  we  adopt  some  notations  of  [CheOl].  The  interrogating  pulses  emitted  by  the  antenna  of  a 
SAR  are  taken  as  linear  upchirps  of  the  form: 

P(t)  =  A(t)eiu°\  where  A{t)  =  Xr{t)eiat2 .  (3.1) 

In  formula  (2.13),  yT(f)  is  the  indicator  function  of  the  interval  of  duration  r: 


Xr(t) 


1,  (6  [— r/2,r/2], 

0,  otherwise, 


and  a  =  B/(2t),  where  B  is  the  bandwidth  of  the  chirp.  Accordingly,  the  instantaneous  frequency  of  the 
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chirp  is  given  by 

Bt 

w(f)=w0H - t  €  [-t/2,t/2],  (3.2) 

T 

where  u> o  is  the  center  carrier  frequency.  The  modulating  function  Apt)  in  formula  (3.1)  is  assumed  slowly 
varying  compared  to  the  fast  carrier  oscillation  elu>ot. 

Suppose  first  that  the  antenna  is  a  motionless  point  source  located  at  a:  £  R3.  Then,  the  propagating 
field  due  to  the  emitted  chirp  (3.1)  is  given  by  the  standard  retarded  potential: 


1  P(t  —  \z  —  x\/c) 
47t  |  z  —  x  | 


(3.3) 


Let  us  assume  that  the  imaged  terrain,  which  is  also  motionless,  is  characterized  by  the  variable  refraction 
index  n  =  n(z).  Under  the  first  Born  approximation  [BW99,  Section  13.1.4],  scattering  is  linearized  so  that 
the  terrain  is  interpreted  as  a  secondary  waves’  source  due  to  the  incident  field  ip(t,  z)  of  (3.3): 


1  -  n2{z)  d2ip  def 

c2  dt2  ~  V[Z)  dt2  ' 

Consequently,  the  scattered  field  at  the  location  x  and  time  t  is  given  by  the  Kirchhoff  integral: 

=  i  IIJ  ^V7|§|((-|i-z|/c.z)<iz.  (3.4) 

\x—z\<.ct 


As  the  amplitude  Apt)  in  (3.1)  is  slowly  varying,  it  can  be  left  out  when  differentiating  the  incident  field 
(3.3)  for  substitution  into  (3.4),  which  yields: 


d2<p 

lit2 


(M) 


P{t-  I Z  -  x\/c) 
47r  \z  —  x\ 


(3.5) 


Consequently, 


fff  _ "M _ p(t 

16tt2  JJJ  \x-z\\z-x\ 

\x—z\<.ct 


\x—z\<.ct 


\x  —  z\\z  —  x\ 


A  (t 


\x  —  z\J  c  —  |  z  —  x\/c)dz 


\X  -z\/c-\z-  x\/c)  e^o{t-\x-z\/c-\z-x\/c)dz 


(3.6) 


Formula  (3.6)  indicates  that  the  scattered  field  ip(t,  x)  is  obtained  by  applying  a  Fourier  integral  operator 
to  the  target  reflectivity  function  v{z),  see,  e.g.,  [CN04,CB08].  The  key  goal  of  radar  imaging  is  to  build  an 
(approximate)  inverse  to  this  operator  and  thus  reconstruct  v(z)  from  the  observed  scattered  field. 

The  approximate  inversion  is  done  in  two  stages;  first,  by  applying  a  matched  filter  to  the  field  ip,  and 
then  by  accumulating  the  information  due  to  multiple  interrogating  pulses  (3.1)  emitted  from  and  received 
at  different  locations  on  the  orbit.  As  one  will  see  from  our  subsequent  analysis  (a  more  detailed  derivation 
can  be  found  in  [Che01,CN04,CB08]),  this  procedure  resembles  the  application  of  an  adjoint  operator,  which 
would  have  coincided  with  the  true  inverse  if  the  mapping  (3.6)  was  a  standard  Fourier  transform  (a  unitary 
operator). 

The  matched  filter  is  defined  as  follows.  Assume  that  there  is  a  point  scatterer  at  a  reference  location 
y,  then  the  resulting  field  at  ( t ,  x)  is  obtained  by  substituting  u{z)  =  S(z  —  y)  into  formula  (3.6): 


U)2  P{t-\x-y\/c-\y-x\/c) 
167 r2  \x  —  y\\y  —  x\ 


(3.7) 


The  filter  is  essentially  a  complex  conjugate  of  ip i  given  by  (3.7);  for  simplicity,  the  constant  factor  — Wg/167r2, 
as  well  the  entire  denominator,  which  is  a  slowly  varying  function  (compared  to  the  fast  oscillation  elulot), 
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are  disregarded.  What  remains  is  merely  P(t  —  \x  —  y\/c  —  \ y  —  x\ / c) ,  where  the  overbar  denotes  complex 
conjugation.  The  application  of  this  filter  yields  the  image: 

I{y)  =  J  P(t- \x  -  y\/c- \y  -  x\/c)ip(t,  x)dt  (3.8) 

=  -  ^2  JJJ  J  P(t~  \x  -  y\/c-\y-  x\/c)P(t-\x  -  z\/c-\z  -  x\/c)dt  _  _  dz, 

\x-z\<cp - - - ' 

W(y,z) 

where  we  have  changed  the  order  of  integration  after  substituting  expression  (3.6)  for  ip{t,  x).  The  interior 
integral  that  we  denote  W ( y ,  z)  in  formula  (3.8)  is  called  the  point  spread  function,5  see  [CheOl].  Up  to  a 
slowly  varying  denominator,  the  point  spread  function  W{y,  z)  yields  the  image  of  a  point  scatterer  located 
at  z,  i.e. ,  it  is  the  field  due  to  a  unit  magnitude  delta-function  at  z  processed  with  the  matched  filter  P(  •  ). 

The  antenna  of  a  real  SAR  sensor  emits  and  receives  a  series  of  chirps  (3.1)  as  the  satellite  travels  along 
the  orbit.  Let  us  therefore  consider  a  sequence  of  emitting  times  and  locations  (tn,  xn)  and  the  corresponding 
sequence  of  receiving  locations  xn.  For  each  n,  we  build  the  point  spread  function  following  (3.8): 

Wn(y,  z)  =  J  P(t-t„-  \xn  -  y\/c—  \y  -  xn\ /c)P(t  -tn-  \xn  -  z\/c-  \z  -  xn\/c)dt.  (3.9) 

The  generalized  ambiguity  function  of  a  SAR  system  takes  into  account  the  information  from  multiple 
interrogating  pulses  by  summing  up  the  corresponding  contributions  (3.9): 

W(y,  z)  =  J2  d(z,  xn)Wn(y,  z).  (3.10) 

n 

The  factor  i 9(z,  xn )  under  the  sum  in  (3.10)  determines  the  range  of  summation.  It  comes  from  the  directivity 
pattern  of  the  antenna,  because  in  real-life  settings  the  antenna  is  never  a  point  monopole  with  isotropic 
radiation.  In  fact,  the  actual  antenna  emits  a  rather  narrow  beam,  see  [CheOl],  and  the  quantity  i9(z,xn) 
can  be  approximated  as  follows: 


..  I  1,  if  the  target  z  is  in  the  beam  emitted  from  xn, 

v(z,  x  )  =  < 

10,  otherwise. 

We  will  later  see  (Section  3.3)  how  to  determine  whether  a  given  target  is  in  the  antenna  beam  or  not.  We 
will  also  see  that  it  is  precisely  the  summation  (3.10)  that  can  help  interpret  the  processing  of  the  received 
field  ip(t,  x)  in  the  sense  of  the  (discrete)  inverse  Fourier  transform  (Sections  3.3  and  3.4). 

The  generalized  ambiguity  function  (3.10)  is  a  key  construct  that  allows  one  to  analyze  the  performance 
of  a  SAR  system.  Indeed,  it  yields  the  image  of  an  ultimately  simplified  target,  which  is  a  point  scatterer, 
whereas  other,  more  complex,  targets  can  be  assembled  of  point  scatterers.  Of  central  importance  in  terms 
of  performance  is  the  resolution  of  a  radar,  i.e.,  its  ability  to  distinguish  between  two  different  targets  located 
at  a  certain  distance  apart.  If  we  could  make  the  generalized  ambiguity  function  (3.10)  equal  to  the  delta- 
function,  W(y,  z)  =  5(y  —  z),  then  the  radar  would  have  had  an  ideal  resolution.  In  reality,  however,  this 
is  never  achieved,  but  the  sharper  (i.e.,  narrower)  the  peak  that  W(y,z)  has  when  the  reference  location  y 
approaches  the  actual  target  location  z,  the  better. 

Note  that  in  most  studies  of  monostatic  synthetic  aperture  radar  resolution,6  the  receiving  locations  xn 
are  taken  the  same  as  the  emitting  locations  xn ,  see,  e.g.,  Section  2.3.1.  This  corresponds  to  the  start-stop 
approximation,  when  the  antenna  emits  the  interrogating  pulse  and  receives  the  scattered  response  while  at 
the  same  position,  and  then  continues  to  move  along  the  orbit  to  the  next  position.  In  this  section  though, 
we  have  intentionally  allowed  for  the  receiving  and  emitting  locations  to  differ  in  the  expressions  for  the 
point  spread  function  (3.9)  and  the  generalized  ambiguity  function  (3.10).  This  will  let  us  take  into  account 
the  actual  antenna  displacement  during  the  pulse  round-trip  time  between  the  orbit  and  the  ground,  and 


*Recall  that  botli  I (y )  and  W(y,  z)  also  depend  on  the  emitting  and  receiving  locations  x  and  x,  respectively. 
(,Bistatic  systems,  unlike  monostatic,  have  separate  transmitting  and  receiving  antennas,  see,  e.g.,  [CGM95,  page  18]. 
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thus  analyze  the  effect  of  the  start-stop  approximation  on  the  quality  of  the  image,  see  Section  3.3.  Note 
also  that  formulae  (3.9)  and  (3.10)  assume  that  both  the  antenna  and  the  target  are  motionless  when  the 
pulse  is  emitted  and  received.  To  account  for  the  motion  and  hence  for  the  Doppler  effect,  the  expressions 
for  the  point  spread  function  and  the  generalized  ambiguity  function  will  be  modified  with  the  help  of  the 
Lorentz  transform,  see  Section  3.4. 

3.3.  Antenna  displacement.  For  the  analysis  of  SAR  resolution,  the  generalized  ambiguity  function 
(3.10)  can  be  conveniently  decomposed  into  the  product  of  two  factors,  transverse  and  longitudinal,  that  will 
be  responsible  for  the  range  and  azimuthal  (i.e. ,  cross-range)  resolution,  respectively,  see  [CheOl].  To  do  so, 
we  first  use  formulae  (3.1)  and  (3.9)  and  recast  expression  (3.10)  as  follows: 


W(y,  z)  =  Yl  ^(z>  x")  / -tn-  \xn  -  y\/c  -| y-  xn\/cy“A\*n-y\/c+\y-*n\/c) 

n  ^ 

xA(t  -tn-  |*n  -  z\/c  -  \z  -  xn\ /c)e-<"°(|s”-*l/c+|z-iB"|/c)dt. 


(3.11) 


Next,  we  take  the  new  integration  variable  as  t  —  tn  in  each  term  of  the  sum  (3.11)  and  denote  it  by  t  again. 
Then,  we  realize  that  neither  A(t  —  \xn  —  y\/c  —  \y  —  xn\/c)  nor  A(t  —  \xn  —  z\/ c—  \z  —  xn\/c )  depends  on 
n  explicitly,  except  for  the  dependence  via  xn  and  xn.  The  latter  is  weak,  because  we  are  assuming  that  the 
distance  from  the  orbit  to  the  ground  is  much  larger  than  the  distance  between  the  successive  emitting  (or 
receiving)  positions  of  the  antenna  on  the  orbit.  Hence,  | y  —  xn\ / c,  etc.,  are  slowly  varying  functions  of  n, 
and  on  top  of  it  we  have  another  slowly  varying  function  A(  •  ).  Therefore,  the  amplitudes  A(  •  )  and  A(  • ) 
can  be  taken  outside  the  sum  (3.11),  and  we  obtain: 


W(y,z)K  (E  $(z?  Xn^eiu0{\x"--y\/c+\y-xn\/c-\xn-z\/c-\z-xn\)  ^ 


WA(y,z) 


(  J  A{i-  \x°  -  y\/c—  I  y  -  x°\ /c)A{t-  \x°  -  z\/c -  \z  -  x°\/c)dt ), 


WR{y,z) 


where  we  have  replaced  the  superscript  “n”  by  “0”  in  the  factor  WR(y,z)  for  definiteness.  Consequently, 
the  generalized  ambiguity  function  can  be  (approximately)  represented  as  a  product  of  the  azimuthal  factor: 


WA{y,z)  =  YJV{z^nYul° 

n 


—  y\/  c+\y—xn  |  /  c—\xn  —  z  \  /  c—\z  —  xn  |/  c) 


(3.12) 


and  the  range  factor: 

WR(y,  z)  =  J  A(t  -  \x°  -  y\/c  -  | y  -  x°\/c)A(t  -  \x°  -  z|/c  -  \z  -  x°\/c)dt.  (3.13) 

We  will  see  that  the  factor  WA{y,z)  of  (3.12)  is  responsible  for  the  resolution  in  the  longitudinal,  or 
azimuthal,  direction  (i.e.,  along  the  orbit),  whereas  the  factor  WR(y,z)  of  (3.13)  is  responsible  for  the 
resolution  in  the  transverse,  or  range,  direction  (perpendicular  to  the  orbit). 

3.3.1.  Azimuthal  resolution.  Hereafter,  we  will  denote  individual  Cartesian  coordinates  by  sub¬ 
scripts:  x  =  (xi,  X2,  X3),  y  =  (2/1, 2/2)  2/3),  z  =  (zi,  Z2,  Z3).  With  no  loss  of  generality  we  can  assume  that  the 
orbit  is  parallel  to  the  Cartesian  direction  “1”  (recall,  we  are  considering  only  short  stretches  of  the  orbit  and 
hence  can  neglect  its  curvature).  To  analyze  the  azimuthal  resolution  we  will  assume  that  both  the  target 
z  and  the  reference  point  y  are  located  on  the  Earth’s  surface  at  the  same  distance  Ro  from  the  orbit,  see 
Figure  3.1.  Let  k  denote  the  angle  between  the  velocity  v  of  the  satellite  (positive  x\  direction)  and  the 
direction  of  the  antenna  beam  (i.e.,  its  central  line).  If  k  <  7r/2,  then  the  antenna  is  pointing  forward;  if 
k  >  7t/2,  then  the  antenna  is  pointing  backward;  and  the  value  k  =  n/2  corresponds  to  broadside  imaging. 

Let  us  denote  the  center  of  the  antenna  beam  footprint  on  the  Earth’s  surface  by  q  =  (<71,  q2,  53).  Clearly, 
we  have:  q\  =  X\  +  i?ocota:,  see  Figure  3.1.  Hereafter,  we  will  assume  with  no  loss  of  generality  that  both 
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Fig.  3.1.  SAR  imaging  schematic  for  the  squinted  stripmap  mode  with  the  antenna  pointing  either  forward  or  backward: 
^forward  7r/2  and  /^backward  >  7r/2. 


the  target  z  and  the  reference  point  y  are  close  to  the  center  of  the  footprint  q  so  that  \qi  —  z±\  <C  Ro 
and  \q-[  —  yi\  <C  Ro-  Then,  we  can  use  the  law  of  cosines  and  the  first  order  Taylor  expansion  of  the  square 
root  to  obtain  approximate  expressions  for  the  distances  that  appear  in  the  exponent  in  formula  (3.12),  see 
Figure  3.1.  For  example,  we  have: 


z  -  x  = 


TD  2  D 

0  +  (Zl  -  <tf)2  -  2-^5_(*1  -  tf)  cos(7r  -  k) 

SU1K 


sin2  k 


Ro 


sin  ac 


/1  (  (zi  —  q™)2  sin2  k  (  0(^i  —  q™)  sin  r  cos  r 

1  +  TTq  + 


Rl 


Ro 


(3.14) 


Ro  ,  (Zl  -  qi)2smn 

- -  +  - -  +  21  -  Ql  COS  K, 

sin  k  2Ra 


and  similarly  for  the  remaining  three  quantities:  |y  —  x"|,  —  y\,  and  \xn  —  z\. 

Expression  (3.14)  and  analogous  expressions  for  other  distances  need  to  be  substituted  into  formula 
(3.12).  Before  doing  that,  let  us  determine  the  range  of  summation  in  (3.12).  Let  L  be  the  size  of  the 
antenna  in  the  direction  parallel  to  the  orbit,  and  let  A  =  27rc/wo  denote  the  carrier  wavelength.  Typically 
we  have  A  <  L,  in  which  case  the  antenna  emits  a  narrow  beam  with  the  angular  width  of  approximately 
2X/L,  see  [CheOl]  and  [Tsy09a,  Appendix  B].  Then,  it  is  easy  to  see  (from  Figure  3.1)  that  the  target  z  is 
in  the  beam  as  long  as 


\zi 


Qi\< 


Ro 
sin  k 


1 

sin  k 


Ro  A 
sin2  k  L 


Consequently,  we  can  write: 

ti(z,xn)  = 


if  Z\  -  <  g"  <  Zl 

L  L  sinz  k  — 

ARp 


ARq 
L  sin2  k,  ’ 


0,  if  qV-  <  Zi  -  T  AK(1  or  q?  >  Zl  +  ,  XR(1 

’  1  L  sinz  k  ±  L  sinz 


(3.15) 


Let  Acci  be  the  distance  along  the  orbit  between  the  successive  emissions  of  pulses  so  that  a:"  =  nAaq.  Let 
us  also  recall  that  q\l  =  a;"  +  Ro  cot  k.  see  Figure  3.1,  and  introduce 


2XRo 

AaqLsin2  k 


and 


N0  = 


Zl  —  Ro  cot  K 
Axi 


N  = 
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where  [  •  ]  stands  for  the  integer  part.  Then,  according  to  (3.15)  we  can  re-write  the  sum  (3.12)  as  follows: 


No  +  N/2 


WA{y,z)=  Y 


oiu0(\xn -y\/c+\y-xn\/c-\xn -z\/c-\z-xn\/c) 


(3.16) 


n=N0-N/2 


Next,  we  substitute  expression  (3.14)  and  the  like  into  the  exponents  in  formula  (3.16).  In  doing  so,  we  take 
the  terms  that  do  not  depend  on  n  outside  the  sum: 

lVo+JV/2 

WA{y,z)  = 

n=N0-N/2 


Nq+N/2 

=  e'“0(o»fa-2,.l)e  (2y  l  -  2z\ ) 

n=N0-N/2 

No+N/2 

=  e  iu°  r  "  Y  _ 

n=N0-N/2 


(3.17) 


We  have  now  arrived  at  the  point  when  we  need  to  identify  a  relationship  between  the  quantities  g"  and  g™ 
or  equivalently,  between  x\  =  gi  —  Ro  cot  k  and  X\  =  g-|  —  Rq  cot  k.  For  each  n,  x™  denotes  the  first  coordinate 
of  the  location  at  which  the  n-th  pulse  is  emitted,  whereas  x"  is  the  first  coordinate  of  the  location  at  which 
the  scattered  signal  is  received.  Clearly,  g"  —  g"  =  x"  —  x”  =  vT,  where  v  =  |v|  is  the  speed  of  the  satellite 
and  T  =  T(z)  is  the  round-trip  travel  time  between  the  antenna  and  the  target.  To  find  T,  let  us  denote  by 
7  the  angle  between  the  velocity  v  and  the  direction  z  —  x  from  the  antenna  x  to  the  target  z\  it  is  easy  to 
see  from  Figure  3.1  that  cot  7  =  ( z\  —  x\)/Rq.  Using  the  law  of  cosines,  we  can  then  write: 


T(z )  =  — - —  +  -\A z  —  x\2  +  ( vT )2  —  2| z  —  x\ (vT)  cosy 

cc 

=  +  -y/\z  -  x\ 2  +  (vT)2  -  2(Zl  -  Xl)(vT). 

c  c 


(3.18) 


This  is  a  quadratic  equation  with  respect  to  T.  Solving  it,  we  obtain: 


rr,^  _  o \z  -  A/c  -  {zi  -  X\) / c{v / c) 
1  yz )  z 

1  —  /  cz 


(3.19) 


Then,  combining  (3.19)  with  (3.14)  and  disregarding  the  v2 /c2  term  in  the  denominator,  we  can  write: 


Qi  ~q?=vT  =  2 


v  Rq  v  (27  —  Qij  sin  k  v 


c  sin  k  c 


Ro 


+  2-(zi  -  g")  cos  k  -  2  —  (27  -  x"). 


(3.20) 


In  the  previous  expression,  we  can  also  drop  the  second  term  on  the  right-hand  side  as  it  is  much  smaller 
than  the  third  term.  Consequently,  we  can  continue  equality  (3.17)  as  follows: 


WA{y,z)=e  °r 


No+N/2 

(2yi-2z1)e^§^L(2yl-2zl)  ^  ^ (Zl-Vl)(2q/ +vT) 

n—No  —  N/2 


No+N/2 

=  e  (2y1-2z1 )  e  -  (2y l -2z2r +2(Zl -Vl )  f  ( ^  +Zl  cos  *))  ^  & '"°  °‘en -  2(Zl -yi  )(1-  g  cos  K)q? 

n—No  —  N/2 


where  we  have  taken  the  factors  that  do  not  depend  on  n  outside  the  sum,  and  besides,  we  have  neglected 
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the  last  term  of  (3.20)  in  the  exponent  as  it  was  only  a  ~  v2/c2  correction  to  the  zeroth  order  term.  Then, 

Wa{V  z )  =  e*"°  l°“  K  (2y1-2z1)  K  (yi-zi+(zi-yi)(^(inr^+zi  cosk)+2(1-^  cos  n)R0  cot  ft)) 

No  +  N/2 

x  ^2 


n=N0-N/2 


-(2y1-2z1) (yi~zi+(zi-?/i)(iKa^+zi  cos  «)+2(l-f  cosk)(_R0  cot  k+IV0  Axx))) 


(3.21) 


N/2 

E 

^=  — 1V/2 


x  ^  g  '^flpc  "2(zi-j/i)(l-l  cosK)riAx1 


Two  observations  can  immediately  be  made.  First,  the  factor  in  front  of  the  sum  in  formula  (3.21)  has 
magnitude  one  and  can  be  left  out  of  consideration  when  analyzing  the  SAR  resolution.  Second,  the  effect  of 
the  antenna  displacement  during  the  round-trip  travel  time  between  the  orbit  and  the  ground  is  accounted 
for  by  the  factor  (1  —  -  cos  ft)  in  the  exponent  under  the  sum  in  formula  (3.21).  To  actually  estimate  the 
azimuthal  resolution,  one  needs  to  compute  the  sum: 


N/2 

WA(v,z)  oc  e^^2izi-Vl)(1-cCOSK)nAxi. 

n=-N/2 


(3.22) 


This  sum  can  be  interpreted  as  the  discrete  (inverse)  Fourier  transform  which,  for  example,  is  encountered 
frequently  in  the  analysis  of  finite-difference  approximations  (see,  e.g.,  [RT07]),  as  well  as  in  a  variety  of 
other  contexts.  Sum  (3.22)  is  a  geometric  sequence  and  can  therefore  be  evaluated  easily.  Again,  up  to  a 
factor  of  magnitude  one  we  have: 


WA{y,z)  oc 


sin 

-  Vl)(l  -  l  cos  k)(N  +  l)AXl 

sin 

’!^(*1-l/1)(l-§cos«)Aa;1' 

(3.23) 


Let  us  further  assume  that 

i  i  ..  Roc  R0X 

zi  -  yi  <  — : — x —  =  x — : — x — > 

ujq  sin  ft AX\  27T  sin  ftAaq 

which  simply  means  that  for  the  purpose  of  conducting  the  analysis  we  would  like  the  reference  point  to  be 
sufficiently  close  to  the  target.  In  this  case  the  argument  of  the  sine  function  in  the  denominator  of  (3.23)  is 
small  and  we  can  write: 


WA(y,  z)  oc  (N  +  l)sinc 


wo  sin  ft 
RqC 


( zi  -  yi)  (l 


-  COS  ft 

c 


{N+l)AXl 


(3.24) 


The  sine  function  (3.24)  is  what  one  typically  gets  instead  of  the  ideal  5(z\  —  y\)  when  analyzing  radar 
performance.  The  sharper  the  maximum  of  the  sine,  i.e. ,  the  closer  it  is  to  the  delta- function,  the  better  the 
resolution  of  the  SAR,  i.e.,  the  minimum  distance  between  the  two  point  targets  that  the  radar  can  still  tell 
apart.  Therefore,  it  is  natural  to  estimate  the  resolution  as  the  width  of  the  main  lobe  of  the  sine: 7 


(Ai  -  2/i )  ( 


1  — -  COS  ft  )  = 


2ttR0c 


RnX 


L  sin  ft 


wq  sinft(7V  +  l)Acci  sin  k(N  +  1) Ax i 


so  that  we  can  write  (v/c  <C  1): 


zi  -  2/i 


L  sin  ft 


^1  +  -  cos  kJ 


(3.25) 


7Sometimes,  it  is  taken  as  half  this  width. 
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If  not  for  the  additional  factor  (l  +  -cos/t),  formula  (3.25)  would  yield  a  standard  expression  for  the  az¬ 
imuthal  SAR  resolution  in  the  framework  of  the  start-stop  approximation  (in  the  squinted  stripmap  mode). 
The  additional  factor  disappears  when  k  =  7t/2.  In  other  words,  in  the  case  of  broadside  imaging,  the 
start-stop  approximation  does  not  affect  the  azimuthal  resolution  of  a  SAR  sensor  at  all.  Otherwise,  when 
cos  k  ^  0,  the  foregoing  analysis  indicates  that  the  actual  resolution  may  be  slightly  worse  than  that  obtained 
with  the  help  of  the  start-stop  approximation.  The  deterioration,  however,  is  minute,  on  the  order  of  v/c, 
which  is  about  2.7 •  10-5  for  low  orbit  satellites,  and  which  can  be  safely  disregarded  for  all  practical  purposes. 

3.3.2.  Range  resolution.  To  analyze  the  range  resolution,  we  substitute  the  envelope  A  of  the  chirp, 
see  (3.1),  into  formula  (3.13): 


where 


and 


WR{y,z) 


^max 


max{(|i°  — 


Train 


min{(|cc0  — 


gia(t—  \x°  —  y\/c—  \y— x°\/c)2  \x°  —  z\/c—  \z— x°\/c)2  (3.26) 

y\/c+\y  -  x°|/c),  (|i°  -  z\/c+\z  -  x°\/c)}  (3.27a) 

y\/c+\y  -  x°\ /c),  (|i°  -  z\/c+\z  -  x°\ /c)}.  (3.27b) 


Formula  (3.26)  holds  provided  that  At  '='  rmax  —  rm;n  <  r,  otherwise  the  integral  is  equal  to  zero.  If  it  is 
not  zero,  then  from  (3.26)  we  find: 

WR(y,z)  =  e1. 


fT/2+T-a, 


'  —t/2  +TItI 


^2iat{(\x0-y\/c+\y-x0\/c)-(\x0-z\/c+\z-x°\/c))  ^ 


(3.28) 


The  factor  in  front  of  the  integral  in  (3.28)  can  be  disregarded  because  in  has  magnitude  one.  The  integral 
itself  needs  to  be  evaluated  for  the  case  when  the  reference  location  y  is  close  to  the  target  z\  then,  the 
length  of  the  interval  of  integration,  r  +  rm ;n  —  rmax  =  r  —  At,  will  also  be  close  to  r.  The  integration  yields: 

WR(y,z)  oc  — 


2ia((\x°  —  y\  +  \y  —  x°\)  —  (\x°  —  z\  +  \z  —  ®°|)) 

t /2+rmin 


x  e 


2ioct((\x°  —  y\/  c-\-\y—x°\/  c)  —  (\x°  —  z\/  c-\-\z—x°\/  c)) 


—  r/2+rm 


or  equivalently  (up  to  another  factor  of  magnitude  one): 


WR{y,z)  oc 


2m((|a:0  —  y\  +  \y  —  *°|)  —  (|*°  —  z\  +  \z  —  a:0|)) 


x  e 


2iat{{\xa -y\/c+\y-x°\/  c)-(\x° -z\/  c+\z-x°\/c)) 


r/2— Ar/2 
-r/2+Ar/2 


2  CT 

B{( \x°  -  y\  +  \y  -  a:0|)  -  (|x°  -  z\  +  \z  -  a:0|)) 


x  sm 


/  ^((l^0  —  y\  +  \y  —  #°|)  —  (\x°  —  z\  +  \z  —  x°|))  BAt2\ 


v 


2c 


± 


2r  )' 


(3.29) 


where  B  is  the  bandwidth  of  the  chirp,  see  (3.2),  and  the  sign  plus  or  minus  under  the  sine  function  in  (3.29) 
depends  on  what  actual  value  delivers  the  maximum  and  minimum  in  formulae  (3.27a)  and  (3.27b).  In  any 
event,  as  we  are  primarily  interested  in  the  case  A r/r  <  1,  we  can  disregard  the  second  term  under  the  sine 
and  obtain: 


WR(y,  z)  (xt  ■  sine 


^R((|^°  -  y\  +  \y_ -  a:0|)  -  (|k°  -  z\  +  \z  -  a;°|))^ 


(3.30) 
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If  the  receiving  location  x°  coincides  with  the  emitting  location  a:0,  which  corresponds  to  the  start-stop  ap¬ 
proximation,  then  formula  (3.30)  yields  a  standard  expression  for  the  range  resolution,  i.e. ,  for  the  resolution 
of  the  distance  between  the  antenna  and  the  target.  It  is  given  by  the  full  width  of  the  main  lobe  of  the  sine, 
which  is  twice  its  half-width,  and  the  latter  is  obtained  by  setting  the  argument  of  the  sine  equal  to  n: 


\x°  y |  -  \x°  -  z\  = 


2l TC 


(3.31) 


If  the  emitting  and  receiving  locations  are  different,  then  instead  of  (3.31)  we  have: 

(l*°  ~y\  +  \y-x°\)~  (\x°  -z\  +  \z-  *0|))  =  (3.32) 

Formula  (3.31)  can  be  interpreted  as  if  the  target  z  and  the  reference  location  y  lie  on  two  concentric  spheres 
centered  at  x°,  and  these  two  spheres  can  be  told  apart  if  their  radii  differ  by  more  than  2irc/ B.  Likewise, 
formula  (3.32)  shall  be  interpreted  as  if  the  target  z  and  the  reference  location  y  lie  on  two  ellipsoids  of 
revolution  with  common  foci  x°  and  x°,  and  these  two  ellipsoids  can  be  told  apart  if  the  corresponding  sums 
of  distances  to  the  foci  (each  sum  is  equal  to  twice  the  major  semi-axis)  differ  by  more  than  47r c/B. 

Hence,  the  eccentricity  of  the  ellipsoids  in  (3.32)  will  provide  a  measure  of  the  difference  between  the 
resolution  obtained  using  the  start-stop  approximation  and  the  resolution  obtained  with  no  use  of  this 
approximation.  It  turns  out  that  this  eccentricity  is  extremely  close  to  one,  i.e.,  that  the  ellipsoids  are 
almost  indistinguishable  from  spheres.  Indeed,  if  the  major  semi-axis  is  estimated  as  a  =  Rq  ~  lOOOfcm, 
then  the  minor  semi-axis  is  b  =  \/ Rq  —  ( vT/2 )2  «  R0(l  —  (vT)2 /8Rq)  ~  Ro(l  —  2- 10-10),  because  vT  ss  40m. 
With  the  chirp  bandwidth  B  >  10 7Hz,  the  resolution  given  by  either  (3.31)  or  (3.32)  appears  to  be  on  the 
order  of  meters,  whereas  the  correction  due  to  the  eccentricity  will  be  fractions  of  a  millimeter.  Clearly, 
it  can  be  disregarded,  which  altogether  means  that  in  practice  the  displacement  of  the  antenna  during  the 
round-trip  travel  time  between  the  orbit  and  the  Earth  does  not  affect  the  resolution  of  a  spaceborne  SAR 
sensor. 


3.4.  The  Doppler  effect.  In  Section  3.3,  we  do  take  into  account  that  the  antenna  changes  its  position 
as  the  signal  it  emits  travels  to  the  target  and  back.  We  still  assume  though  that  the  antenna  is  motionless 
during  the  time  the  signal  is  emitted  and  the  scattered  response  received.  In  reality,  however,  the  antenna 
is  moving  continuously.  Therefore,  the  expressions  for  the  emitted  and  scattered  field  that  we  employed 
for  constructing  the  generalized  ambiguity  function  in  Section  3.2  may  themselves  need  to  be  corrected.  In 
particular,  expression  (3.3)  yields  the  Held  emitted  by  a  motionless  antenna;  it  will  need  to  be  modified  to 
account  for  the  actual  motion.  Likewise,  expression  (3.7)  yields  the  scattered  field  in  the  original  motionless 
coordinate  system;  it  will  need  to  be  modified  because  the  field  is  received  by  a  moving  antenna. 

In  this  section,  we  will  use  the  Lorentz  transform  to  address  the  fact  that  the  antenna  is  moving.  We  will 
see,  in  particular,  that  once  the  motion  of  the  wave’s  emitter  and/or  receiver  is  properly  accounted  for,  the 
observable  frequency  of  the  signal  changes.  This  phenomenon  is  known  as  the  Doppler  frequency  shift,  see, 
e.g.,  [LL75]  for  electromagnetic  waves  and  [LL87]  for  acoustics.  In  the  literature  on  SAR,  this  frequency  shift 
is  sometimes  used  to  explain  the  mechanism  of  azimuthal  resolution,  because  the  sequence  of  signals  that 
illuminate  the  target  from  different  emitting  positions  along  the  orbit  can  be  interpreted  as  a  linear  chirp, 
see,  e.g.,  [FL99].  A  relation  between  the  Doppler  approach  and  that  often  used  in  tomography  has  also  been 
established  [MOJ83,  AM89].  However,  in  the  rigorous  analysis  based  on  the  generalized  ambiguity  function, 
the  Doppler  effect  is  often  disregarded.  We,  on  the  other  hand,  will  show  that  the  Doppler  frequency  shift 
should  be  included  into  the  definition  of  a  matched  filter  so  that  to  minimize  the  distortions  of  the  image. 

Let  us  assume  that  the  antenna  is  engaged  in  a  straightforward  uniform  motion  along  the  coordinate 
“1”  of  a  Cartesian  coordinate  system,  and  that  the  speed  of  motion  is  v,  see  Figure  3.1.  Of  course,  strictly 
speaking  this  is  not  true  for  a  satellite  orbiting  the  Earth,  but  it  is  a  sufficiently  good  approximation  for 
short  stretches  of  the  orbit.  We  will  denote  the  position  of  the  antenna  by  x  =  x(t)  =  (vt,H,  0),  where  H 
is  the  orbit  altitude,  see  Figure  3.1,  and  we  will  also  assume  that  the  pointwise  antenna  emits  a  signal  P(t). 
To  find  the  corresponding  field  ip(t,  z)  as  a  function  of  the  time  t  and  the  observer’s  (target)  location  z ,  we 
first  change  the  coordinates  so  that  in  the  new  frame  the  antenna  becomes  motionless.  Rising  the  standard 
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notation  for  the  relativistic  square  root:  P  =  —  v2 /c2,  where  c  is  the  propagation  speed,  we  write: 

(3.33a) 
(3.33b) 


t  v/c  Zi 


,  v/c 

c'=  -Tct 


z  1 


The  remaining  two  spatial  coordinates,  z2  and  23,  do  not  change,  so  that  the  overall  transformation  is 
(t,  z±,  Z2,  z3)  1  y  (0,  Cl,  z2,  23).  It  is  known  as  the  Lorentz  transform,  and  its  distinctive  feature  is  that  it 
preserves  the  form  of  the  d’Alembert  operator  [Gar07,  Chapter  3,  Section  3].  In  particular,  the  equation 


1  d  V  d2p  d2cp  d2ip 

c2  dt2  dzf 


dz2 


dz | 


(3.34a) 


=  P(t)S(z  -  x(t))  =  P(t)S(z1  -  vt)S(z2  -  H)S(z3), 


becomes 


1  d2ip  d2p 
c2  dtf2”^ cf  ” 


d2p  d2ip 


dz2 

v/c  Cl 
P  c 


dz2 


)S(fiCl)S(Z2  -  H)8{z3)  =f  mCl,Z2,Z3). 


(3.34b) 


Note  that  if  the  standard  Galileo  transform  was  used  instead  of  (3.33),  then  the  differential  operator  in  (3.34b) 
would  not  have  been  the  same  as  in  (3.34a),  even  though  the  antenna  would  still  have  been  motionless  in 
the  new  frame.  Solution  of  equation  (3.34b)  subject  to  zero  initial  conditions  can  be  easily  found  with  the 
help  of  the  Kirchhoff  integral  in  the  new  coordinates: 


ip{0Xi,z2,z3)  =  — 


f(9  -  p'/c ,  =  1  ~  P/c)/P ) 


p'<.cO 


p' 


4-7T/3 


(3.35) 


where  p'  =  \J (Ci  —  Ci)2  +  ( zi  ~  z'2 )2  +  (23  —  z'3 )2  and  p  =  a/Ci  +  (z2  —  H )2  +  z2.  Note  that  up  to  the 
relativistic  factors  /3,  solution  (3.35)  is  a  standard  retarded  potential  [compare  with  formula  (3.3)]  from  the 
point  source  located  at  (0,  H ,  0)  in  the  new  coordinates,  which  is  expected.  Note  also  that  as  an  alternative, 
solution  (3.35)  could  have  been  obtained  in  a  “brute  force”  way,  by  computing  the  convolution  with  the 
fundamental  solution  of  (3.34a)  in  the  original  coordinates,  see  Appendix  D. 

To  understand  the  Doppler  effect,  solution  (3.35)  needs  to  be  cast  back  to  the  original  coordinates.  Let 
us  denote  r  =  \j z\p  (z2  —  H )2  +  z\.  Then,  according  to  (3.33b)  we  can  write: 


Ci  =  7^(21  -  vt)2  =  —(zl~  2zi vt  +  (vt)2)  =z\  + 


1-/32 


~i  +  m(-2z  ivt  +  (vt)2) 


=  z2  —  2z\vt  + 


1  - 


/ ? 2 


-(zi  —  vt)2  +  (vt)2  =  z2  —  2rvtcosj  + 


P 
l-/?2 
P2 


(zi  -  vt)2  +  (vt)2 


where  as  before  (in  Section  3.3),  7  denotes  the  angle  between  the  velocity  v  (positive  first  coordinate)  and 
the  direction  from  the  antenna  to  the  target,  see  Figure  3.1.  Consequently, 


p2  =  r2  —  2rvtcosj  +  (vt)2  + 


l-/?2 

P2 


(zi  -  Vt)2. 


(3.36) 


The  last  term  in  equality  (3.36)  is  ~  i>2/c2,  where  v/c  itself  is  small,  v/c  <C  1.  The  second  to  last  term  is 
also  quadratic  with  respect  to  v/c  for  the  times  of  interest  t  ~  r/c.  Hence,  for  the  values  of  7  away  from  7r/2 
(so  that  cos  7  0)  these  two  terms  can  be  dropped,  which  yields  the  following  approximate  expression: 


/—z -  I  2vt 

p  k,  \J r1  —  2 rvt  cos  7  =  r  y  1 - cos  7  zz  r  —  vt  cos  7. 


(3.37) 
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Then,  substituting  the  value  of  p  given  by  (3.37)  into  the  solution  (3.35),  using  the  definition  of  the  Lorentz 
transform  (3.33),  and  approximating:  1//3  =  (1  —  v2/c2)-1/2  «  1  +  v2/2 c2  «  1,  we  arrive  at 


1  ^((g(*-^)-|  +  ^CQS7) 
47T/3  p 


1  ^((*-^)  (1  +  yCQS7)) 
47r  p 


(3.38) 


If  the  function  P(t)  is  a  fast  oscillation  with  the  frequency  w0,  then  it  is  also  possible  to  replace  p  by  r  in 
the  denominator  of  (3.38)  because  it  only  amounts  to  a  slow  and  small  (~  v/c)  variation  of  the  amplitude: 


p{t,z) 


1  ^  ((*-f)( l  +  gcos7)) 

47t  r 


(3.39) 


Formula  (3.39)  is  almost  the  standard  retarded  potential,  except  for  the  correction  factor  1  +  -  cosy  inside 
P(  ■  ).  The  presence  of  this  factor  means  that  if  the  field  emitted  by  a  moving  antenna  is  received  at  a 
motionless  target,  then  the  frequency  measured  at  the  target  changes,  instead  of  the  original  wo  it  becomes 
wo  (l  +  -  cosy).  This  implies  the  increase  of  the  frequency  if  the  antenna  is  moving  toward  the  target, 
0  <  y  <  7t/2,  and  the  decrease  of  the  frequency  if  the  antenna  is  moving  away  from  the  target,  tt/2  <  y  <  n. 
The  correction  to  the  frequency  is  linear  with  respect  to  v/c  (since  we  have  dropped  the  quadratic  terms 
~  v2  / c2  in  our  analysis)  and  accordingly,  the  foregoing  phenomenon  is  called  the  linear  Doppler  effect  or 
linear  Doppler  frequency  shift.  Note  also  that  in  the  same  linear  framework  we  can  recast  formula  (3.39)  as 


47 rr 


— hii_y 

1-  fcoe 7/ 


(3.40) 


Expression  (3.40)  can  be  obtained  independently,  after  equation  (3.34a)  has  been  solved  directly  by  comput¬ 
ing  the  convolution  with  the  fundamental  solution,  see  Appendix  D. 

For  the  sake  of  completeness,  let  us  also  comment  briefly  on  the  quadratic  Doppler  effect.  The  last  two 
terms  that  we  have  disregarded  on  the  right  hand  side  of  (3.36)  become  important,  i.e.,  comparable  to  the 
second  term,  if  cosy  ~  v/c  or  equivalently,  7r/2  —  u/c<y<7r/2  +  v/c,  where  v/c  <C  1.  This  means  that 
the  quadratic  effect  can  manifest  itself  only  for  the  broadside  imaging,  when  the  direction  of  the  beam  is 
normal  to  the  orbit.  On  top  of  that,  it  may  only  be  noticeable  in  a  very  narrow  central  part  of  the  beam, 
because  in  practice  v/c  ~  10”5,  whereas  the  angular  semi- width  of  the  beam  X/L  could  be  about  1/30  for 
the  wavelength  A  =  30cm  (corresponds  to  the  carrier  frequency  wo  =  1  GHz)  and  the  antenna  size  L  =  10 m. 
Hereafter,  the  quadratic  Doppler  effect  ~  v2 /c2  ~  10” 10  will  be  neglected. 

To  account  for  the  linear  Doppler  effect  when  building  the  generalized  ambiguity  function,  we  replace 
formula  (3.3)  with  formula  (3.39).  Then,  we  differentiate  c pit,  z)  and  thus  derive  the  sources  of  the  scattered 
held  in  the  framework  of  the  first  Born  approximation  [cf.  formula  (3.5)]: 

3V,  *  ,  W2  (1  +  l  cosy)2  P((i-i)(l  +  |  cosy)) 

dt*{'t,Z)~  r  '  [  ’ 

For  the  point  target  located  at  z,  the  scattered  wave  at  (t,x)  is  a  retarded  potential  centered  at  z  and 
modulated  by  the  function  (3.41)  as  a  function  of  time  [cf.  formula  (3.7)]: 


wg(l+fcosy)2p((t 

167T2 


(3.42) 


When  the  scattered  pulse  (3.42)  travels  back  from  the  target  to  the  antenna,  it  generates  another  contribution 
to  the  Doppler  frequency  shift.  The  corresponding  analysis  can  be  conducted  as  before,  by  going  to  the  frame 
of  reference  in  which  the  receiver  (antenna)  is  motionless.  To  do  so,  we  employ  essentially  the  same  Lorentz 
transform  except  that  as  we  want  the  solution  in  the  form  of  a  retarded  potential  centered  around  the  origin, 
we  need  to  have  the  origin  of  the  new  system  initially  located  at  the  target  z,  rather  than  at  a:(0)  =  (0,  H,  0). 
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As  the  scattered  field  is  received  at  the  antenna  x,  we  will  use  the  variables  x  =  (x\,X2,  £3)  to  write  the 
transform,  whereas  the  location  of  the  target  z  =  (zi,z2,Z3)  will  be  considered  fixed  [cf.  formulae  (3.33)]: 


9 

£1 


t 

0 


v/c  X\  —  ZI 


v/c 


ct 


(3.43a) 

(3.43b) 


Similarly  to  (3.33),  the  remaining  two  spatial  coordinates  do  not  get  changed  so  that  altogether  we  have 
the  transformation  (t,x  1,2:2,  £3)  H >  (9,  £1,  X2,  £3).  Relation  (3.43b)  implies,  in  particular,  that  in  the  new 
coordinates  the  antenna  is  not  moving,  because  x\  =  vt  corresponds  to  £1  =  —Z\/j3  =  const. 

Let  us  introduce  p  =  y/£f  +  (x2  -  z2)2  +  (X3  -  Z3)2  and  r  =  \J (ay  -  Z\)2  +  (x2  -  z2)2  +  (x3  -  z3 )2  = 
\x  —  z\.  Then,  x\  —  z\  =  — r cosy,  see  Figure  3.1,  and  repeating  the  same  derivation  as  led  us  to  formula 
(3.37),  we  arrive  at  the  following  linearized  relation  between  the  distances: 


p  ss  f  +  vt  cos  7. 


(3.44) 


Consequently,  letting  /3  ~  1  and  using  (3.43b),  we  can  write: 

p cosy  «  r  cosy  +  vt  cos2  y  =  —  {x\  —  zi)  +  vt  cos2  y 
«  —  £1  —  vt  +  vt  cos2  y  =  —  £  —  vt  sin2  y 


so  that 


£1  =  —p cosy  —  vt  sin2  y. 


(3.45) 


Next,  using  formula  (3.43a)  and  again  neglecting  quadratic  terms  with  respect  to  v/c,  we  obtain  from  (3.44): 


v  ^ 

f  «  p  —  vt  cosy  =  p  —  vj39  cos  y - o(xi  ~~  zi)  ~  P  ~  v9  cosy.  (3.46) 

cz 

Finally,  to  implement  the  transformation  in  formula  (3.42),  we  have  to  transform  all  times  and  distances  on 
its  right-hand  side  to  new  coordinates.  We  note  that  formula  (3.42)  is  written  so  that  x  can  be  any  location 
in  space.  We,  however,  are  interested  only  in  the  case  when  x  is  the  location  of  the  antenna.  Since  the 
displacement  of  the  antenna  during  the  round-trip  travel  time  of  the  signal  does  not  affect  the  performance 
of  the  radar  (Section  3.3)  we  take:  x  =  a:(0)  =  (0,77,0).  In  this  case,  r  =  \x  —  z\  =  r  and  instead  of  (3.42) 
we  can  write: 


i>(t,x) 


^0  (l  +  c  cosy)2  P  ((t  —  (l  +  f  cosy)) 

16tt2  r'z 


(3.47) 


Substituting  f  given  by  (3.46)  into  formula  (3.47),  and  also  substituting  the  expression  for  t  derived  by 
inverting  (3.43)  and  using  (3.45): 


9  u/c£  1 

t  =  -  + 

0  0  c 


9 


v  p 

—  cosy, 
c  c 


we  have: 


1,X2,X3) 


^(l  +  gcosy )2-P((^¥)  (l+7^7)2) 
167t2  r2 

tag  (1  +  2^  cosy)  P  ((<?  ~  f )  (l  +  2|  cosy)) 
16tt2  r2 

^2  P((o-2j)0  +  2^ cosy)) 

167T2  p2 


(3.48) 
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where  the  last  equality  in  (3.48)  was  obtained  as  before,  by  disregarding  a  slow  and  small  (~  v/c)  variation 
in  the  amplitude.  Formula  (3.48)  indicates,  in  particular,  that  in  the  frame  of  reference  in  which  the  satellite 
is  motionless,  the  observed  frequency  of  the  signal  emitted  by  the  antenna,  scattered  off  the  target,  and  then 
again  received  by  the  antenna  is  lvq  (l  +  2-  cosy)  compared  to  the  original  frequency  ujq.  Otherwise,  if  not 
for  the  coordinate  change  and  the  additional  amplitude  factor,  the  solution  (3.48)  looks  exactly  the  same  as 
it  would  have  looked  if  both  the  antenna  and  the  target  were  motionless  [cf.  formula  (3.7)];  with  —  being 
the  round-trip  travel  time  between  the  antenna  and  the  target.  Hence,  this  scattered  solution  can  be  used 
for  building  the  generalized  ambiguity  function  in  the  new  coordinates. 

3.4.1.  Azimuthal  resolution.  Hereafter,  we  will  follow  the  pathway  of  Section  3.2  for  constructing 
the  generalized  ambiguity  function.  In  doing  so,  we  will  use  solution  (3.48)  instead  of  (3.7)  as  the  central 
“building  block.”  Following  the  analysis  of  Section  3.3,  we  will  also  disregard  the  displacement  of  the  antenna 
from  the  very  beginning.  For  the  azimuthal  factor  we  can  then  write  [cf.  formula  (3.12)]: 

WA{y,  z)=^2  #(z>  xn)e2iu,0^xn~y^1+^  COS7»  Hxn-zl(i+^f -7”))/ci  (3.49) 


where  the  distances  \xn  —  y\  and  \xn  —  z\,  even  though  denoted  the  same  way  as  before,  shall  be  measured  in 
the  Lorentz  transformed  coordinates,  in  which  the  antenna  is  not  moving,8  and  the  notation  for  the  angles, 
7 ™  and  7",  indicates  that  they  depend  on  the  position  of  the  antenna. 

We  are  particularly  interested  in  analyzing  that  part  of  the  exponent  in  formula  (3.49)  which  has  the 
factor  v/c  in  front  of  it: 


2iwo  2v 
c  c 


iujQ  2v 


(h 


x' 


X 


2/|  cosy™ 


|  cosy™) 


Note  that  the  previous  formula  was  derived  with  the  help  of  (3.14).  To  simplify  subsequent  analysis,  we  will 
assume  that  the  imaging  is  broadside,  i.e.,  k  =  7t/2,  see  Figure  3.1.  Then, 

2iuio  2v 


{\xn  -  y\ cosy™  -  |x™  -  z\ cosy™) 


2muo  2v 


i?o(cosy”  -  cosy™)  +  ^ cosy™ 


,n\2 


(Zl  ~  g?) 
2i?o 


■  cos  yz 


Besides,  as  the  antenna  beam  is  narrow,  we  can  also  say  that  both  7™  and  7™  are  close  to  tt/2  so  that 


n  ,  n  Vl  ~  9l  -  n  ,  n  *1  ~  Ql 

cos  7  «  cot  7  =  — — -  and  cos  yz  «  cot  yz  =  — — - 

i?o  Rq 


Hence,  we  have: 


2muo  2v 


c  c 
2iuJo  2v 
RqC  c 


( \xn  -  y\ cosy™  -  |sc™  -  z\ cosy”) 

no* -«)+<»-«?  (*■-*)■ 


2Ro 


2i?n 


(3.50) 


The  key  consideration  regarding  the  previous  formula  is  that  the  fist  term  in  the  brackets  on  its  right-hand 
side  does  not  depend  on  n.  Consequently,  when  expression  (3.50)  is  substituted  back  into  the  exponent  in 
formula  (3.49),  the  corresponding  factor  can  be  taken  outside  the  sum  and  disregarded  because  its  magnitude 
will  be  equal  to  one.  Therefore,  formula  (3.49)  will  eventually  simplify  to  [cf.  formula  (3.17)]: 


No+N/2 

WA{y,z)oc  ^2  e 

ti—Nq  —  N/2 


fl0. 


((zi-yi)qr+g  (yi rq1  ) 


R0 


(3.51) 


8 This  applies  to  all  other  distances  that  appear  in  subsequent  derivations  in  this  section. 
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It  is  clear  that  the  second  and  third  terms  in  the  exponent  on  the  right-hand  side  of  formula  (3.51)  are  much 
smaller  than  the  first  term,  because  in  addition  to  the  small  factor  v/c  they  contain  another  small  factor 
(; yi  —  qi)/Ro  ~  (zi  —  qi)/Ro  ~  A /L.  Dropping  these  two  terms,  we  obtain  [cf.  formula  (3.22)]: 

N 72  2. 

WA(y,z)cx  J2  eT^{zi~Vl)nAxi.  (3.52) 

n=~N/2 


Formula  (3.52)  provides  a  conventional  expression  for  the  azimuthal  component  of  the  generalized  ambiguity 
function,9  see,  e.g.,  [CheOl].  Once  the  sum  on  the  right-hand  side  is  evaluated,  it  yields  the  azimuthal 
resolution: 


zi  ~  Vi  « 


L 

2' 


(3.53) 


which  is  standard  in  the  analysis  of  SAR  performance  (note  that  formula  (3.25)  reduces  to  (3.53)  for  k  =  7t/2). 
Formula  (3.53)  shows  that  if  the  reference  distance  \xn  —  y\  in  the  exponent  in  formula  (3.49)  is  taken  with 
the  Doppler  correction,  i.e.,  if  the  frequency  shift  is  included  in  the  definition  of  a  matched  filter,  then  the 
azimuthal  resolution  of  the  radar  stays  basically  unaffected. 

Let  us  now  see  what  will  change  if  the  matched  filter  is  taken  with  no  Doppler  correction.  In  this  case, 
instead  of  formula  (3.49)  we  have: 


WA{y,  x")e2M\*n-y\-\xn-z\(i+^coSj:))/c 


(3.54) 


and  consequently,  the  additional  terms  in  the  exponent  of  (3.54)  are  [cf.  formula  (3.50)]: 


2 iujQ  2v 


\x  —  z  |  cos 7Z 


2iwo  2u 
RqC  c 


Ro{zi  ~  Qi)  + 


(Zi  -  Qi  )3 
2R0 


(3.55) 


The  first  obvious  difference  that  we  observe  between  formulae  (3.50)  and  (3.55)  is  that  the  leading  term  in 
the  brackets  on  the  right-hand  side  of  (3.55)  depends  on  n.  Hence,  if  we  employ  the  same  argument  as  before 
and  disregard  the  second  term  on  the  right-hand  side  of  (3.55),  then,  instead  of  (3.52),  we  will  arrive  at 


N/2 

WA(y,z)cx  Y.  e^^~y^Ro>Ax' 

n=-N/2  (3.56) 

R^(N +  l)AXl  . 

Formula  (3.56)  implies  that  the  azimuthal  resolution  of  the  radar  per  se  remains  unchanged.  However,  the 
entire  imaged  scene  gets  shifted  along  the  track  (orbit)  by  —  yr-Ro-  For  the  typical  parameters  involved,  the 
magnitude  of  this  spatial  shift  could  be  about  50m. 

3.4.2.  Range  resolution.  For  the  range  part  of  the  generalized  ambiguity  function  we  can  write 
instead  of  (3.26): 


(N  +  l)sinc 


u0  { ,  s  .  2v 


WR{y,z) 


r/2+Tm 


'  —r /2+t 


x  eiw0t(ir  COS7*  — ^  C0B7V)^ 


2 


(3.57) 


where  the  last  factor  under  the  integral  on  the  right-hand  side  of  (3.57)  appears  because  the  Doppler  correc¬ 
tion  shall  be  substituted  not  only  into  the  slowly  varying  envelope  A(t)  but  also  into  the  fast  carrier  oscillation 
e*“° 4  in  the  definition  of  the  chirp  (3.1).  Next,  we  assume  with  no  loss  of  generality  that  yz  =  =  7,  which 


9Formula  (3.22)  yields  a  somewhat  more  complicated  expression  that  accounts  for  the  displacement  of  the  antenna. 
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merely  implies  that  the  reference  location  y  and  the  target  z  are  on  the  same  line  of  sight  from  the  antenna. 
Then,  formula  (3.57)  transforms  into 


WR(y,z) 


(3.58) 


where  for  the  limits  of  integration  we  have  [cf.  formulae  (3.27a),  (3.27b)]: 


2v 

Tmax  =  2  (  1  +  —  cosy  )  max{|a;0  -  y\/c,  \x°  -  z |/c} 


and 


2v 


Tmin  =2  1  +  —  cos 7  min{|xu  -  y\/c,  \xL 

'  c  ' 


z|/c}. 


In  the  exponent  in  formula  (3.58),  we  drop  the  terms  proportional  to  v2 /c2  and  also  disregard  the  terms  that 
do  not  depend  on  t  as  the  latter  only  yield  an  inessential  factor  of  magnitude  one  in  front  of  the  integral: 


WR(y,z)  oc 


Ct/2+t„ 


•  — r /2+r 


*  4( |xu  —  y |  —  \  x'J  —z\) 


(l+^co»7)dti 


(3.60) 


Then,  following  the  same  argument  as  in  the  end  of  Section  3.3,  we  obtain  [cf.  formula  (3.30)]: 


V 


WR(y,  z)  oc  rsinc 

Formula  (3.61)  yields  the  range  resolution: 


f  B{ \x°  -  y\  -  |a:0  -  z\) 


4v 

1  H - cos  7 
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\x°  -  y\  -  \x°  -  Z I  = 


2-7TC 

IT 


,  4v 

1 - cos  7 

c 


(3.61) 


(3.62) 


which  differs  from  the  standard  expression  (3.31)  only  by  the  factor  (l  —  —  cosy),  which  is  very  close  to  one 
(recall,  v/c  ~  10-5).  Moreover,  in  the  case  of  broadside  imaging  we  can  take  7  =  n/2  so  that  (3.62)  reduces 
to  (3.31)  and  the  range  resolution  stays  completely  unaffected  by  the  Doppler  frequency  shift. 

Next,  we  need  to  analyze  the  case  when  the  matched  filter  it  taken  with  no  Doppler  correction.  Then, 
instead  of  (3.60)  we  obtain: 


WR(y,z)  oc 


4(l»° 


J-£|) 


/r  cos  7  j  +iat2  ^  cos  7 


dt. 


(3.63) 


First,  we  notice  that  if  7  =  7t/2  (broadside  imaging),  then  integral  (3.63)  reduces  to  standard  and  yields  a 
conventional  expression  for  the  range  resolution:  \x°  —  y |  —  \x°  —  z\  =  see  formula  (3.31).  Otherwise, 
there  is  an  additional  term  in  the  part  of  the  exponent  in  formula  (3.63)  that  is  linear  with  respect  to  t.  It 
will  cause  an  overall  shift  of  the  entire  imaged  scene  by  \x°  —  z\  —  cosy.  Generally  speaking,  this  value  may 
be  of  the  same  order  of  magnitude  as  that  obtained  for  along-the-track  shift,  see  formula  (3.56);  however, 
for  7  close  to  n/2  (near  broadside  imaging)  it  becomes  small. 

Moreover,  for  7  7^  7r/2  there  remains  a  small  quadratic  term  with  respect  to  t  in  the  exponent  under  the 
integral  on  the  right-hand  side  of  formula  (3.63).  Integrals  of  this  type  were  studied  in  [Tsy09a]  in  the  context 
of  SAR  imaging  through  the  Earth’s  ionosphere.  It  was  shown  that  this  quadratic  term  does  not  affect  the 
resolution  per  se  in  the  sense  that  the  distance  between  the  central  maximum  and  the  first  minimum  of  the 
resulting  ambiguity  function  appears  exactly  the  same  as  half-width  of  the  corresponding  sine.  However, 
the  sharpness  (contrast)  of  the  image  may  be  affected,  because  the  value  of  the  first  minimum  is  not  zero. 
Based  on  the  analysis  presented  in  [Tsy09a],  it  is  possible  to  show  that  for  the  range  factor  given  by  (3.63) 
we  have: 


min  \  WR(y,  z)\ 
max  \WR{y,  z)\ 


ar2  4u 
27t2  c 


cosy  = 


Bt  4v 

tt2  f 


(3.64) 


cos  7. 
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Formula  (3.64)  indicates  that  generally  speaking,  the  degree  to  which  disregarding  the  Doppler  frequency 
shift  may  adversely  affect  the  image  depends  on  the  waveform  used  for  imaging.  More  precisely,  to  reduce 
the  deterioration  of  the  image  sharpness,  one  has  to  minimize  the  product  of  the  pulse  bandwidth  B  times 
its  length  r.  Most  modern  SAR  systems  use  high  range-resolution  pulses  that  have  sufficiently  large  B  [i.e., 
large  chirp  rate  a ,  see  formulae  (3.1)  and  (3.2)];  in  this  case,  the  chirp  duration  must  be  very  short.  It 
turns  out  that  for  the  typical  values  of  the  parameters:  B  =  10 MHz  and  r  =  5  •  10_5sec,  the  degradation 
of  the  image  is  indeed  negligible  because  the  right-hand  side  of  (3.64)  is  ~  10-3cos7-  On  the  other  hand, 
the  so-called  Doppler  SAR  imaging  exploits  the  high  Doppler-resolution  wavetrains  that  approximate  the 
continuous  wave  elLUot  and  are  characterized  by  large  r;  in  this  case,  the  bandwidth  B  must  be  very  narrow. 
Note  that  imaging  using  various  waveforms  has  been  thoroughly  studied  in  [CB08] . 

Let  us  also  note  that  if  the  cubic  terms  were  retained  when  expressions  (3.50)  and  (3.55)  are  substituted 
into  the  exponent  in  formula  (3.49),  we  would  have  arrived  at  a  similarly  small  degradation  of  the  sharpness 
in  the  azimuthal  direction.  We  would  have  also  obtained  a  factor  of  0(1  +  -)  change  in  the  azimuthal 
resolution  itself  that  could,  of  course,  be  disregarded. 

3.5.  Summary.  In  this  part  of  the  report,  we  have  looked  into  whether  or  not  the  start-stop  approxi¬ 
mation  used  for  processing  radar  data  may  affect  the  performance  of  a  SAR  sensor  in  the  case  of  spaceborne 
imaging.  Specifically,  we  have  analyzed  two  effects  that  are  normally  left  out  of  consideration  when  building 
the  generalized  ambiguity  function  of  a  radar.  They  are  the  displacement  of  the  antenna  during  the  pulse 
round-trip  travel  time  between  the  orbit  and  the  Earth’s  surface,  and  the  Doppler  frequency  shift.  Our 
analysis  shows  that  the  displacement  of  the  antenna,  even  though  quite  noticeable  itself,  in  practice  causes 
no  deterioration  of  the  image  resolution.  As  for  the  Doppler  effect,  if  the  frequency  shift  is  included  in  the 
definition  of  a  matched  filter,  then  the  performance  of  the  radar  also  remains  largely  unaffected.  Otherwise, 
the  imaged  scene  gets  shifted,  and  there  may  be  a  small  deterioration  of  the  image  sharpness.  It  should  also 
be  mentioned  that  the  authors  of  [CB08]  refer  to  the  matched  filter  modified  by  the  Doppler  frequency  shift 
as  to  a  phase-corrected  matched  filter. 

A  useful  future  extension  of  the  foregoing  analysis  could  include  taking  into  account  the  curvature  of 
the  orbit. 

Appendix  A.  Propagation  of  radar  pulses  in  the  homogeneous  ionosphere.  We  are  interested 
in  obtaining  a  spherically  symmetric  solution  similar  to  the  retarded  potential  (2.15),  but  for  the  case  of 
dispersive  propagation  governed  by  the  Klein-Gordon  equation  (2.6).  Let  r  =  \z  —  x\  denote  the  radial 
coordinate  in  the  spherical  system  centered  at  s  €  R3,  which  is  the  location  of  the  antenna  (we  are  using 
the  notations  of  Section  2.3.1),  and  let  ip  =  p(t,r)  be  a  spherically  symmetric  solution  of  equation  (2.6). 
Introduce  a  new  function  </>  =  <p(t,r),  such  that  < p(t,r )  =  (f>(t,r)/r,  then  satisfies  the  one-dimensional 

Klein-Gordon  equation: 


d2<j> 

W 


2  d2(j)  2 

C  Qr2  +  Wpe^  —  0’ 


r  >  0. 


(A.l) 


Consequently,  we  need  to  look  into  the  propagation  of  pulses  governed  by  equation  (A.l). 

Assume  that  a  pulse  of  shape  P(t)  is  given  at  r  =  0  (location  of  the  antenna);  for  example,  it  can  be  the 
high-range  resolution  upchirp  (2.13): 


P(t)  =  -^xT(t)ei("“+&)t,  (A. 2) 

47T 

where  wo  is  the  center  carrier  frequency,  B  is  the  bandwidth,  and  r  is  the  duration  of  the  chirp.  The  factor 
1/47T  in  (A. 2)  accounts  for  the  difference  between  the  one-dimensional  and  three-dimensional  delta-functions 
S(r)  and  6(x  1,2:2,  £3)  that  excite  the  pulse  at  the  origin. 

We  will  Fourier  transform  the  pulse  (A. 2)  in  time  and  study  the  propagation  of  individual  frequencies. 
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The  transformation  yields: 


r/2 

t/2 

/  J_e*(wo  +  #F 

y-^dt  -2—  j 

J  47T 

87T2  J 

t/2 

-r/2 
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L<  ( ~  aT)' 
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8-7T2  2yfa 


—  erf 


/  +  cut) 

V  2y/a 


where  a  =  41  is  the  chirp  rate,  and  fj  =  w0  —  oj.  To  simplify  expression  (A.3),  we  need  to  analyze  the 
arguments  of  the  erf  functions.  Denote 


V±  = 


/3  ±  ar  wo  —  w  ±  B/2 


2^ 


(A.4) 


and  consider  several  cases.  First,  let  w  G  [wo  —  B/2,  w>o  +  B/2]  and  in  addition  suppose  that  w  is  not  too  close 
to  either  endpoint,  ojq  —  B/2  or  ujq  +  B/2.  In  other  words,  suppose  that  the  numerator  on  the  right-hand  side 
of  (A.4)  is  of  the  same  order  of  magnitude  as  the  bandwidth,  i.e.,  |wo  —  w±B/2\  ~  B.  Then,  |?j±|  ~  \JBt  / 2, 
and  for  the  typical  values  of  the  parameters  involved  (see  Section  2.3)  this  quantity  is  approximately  15.8. 
In  general,  we  can  say  that  the  value  of  \JBt / 2  is  between  10  to  20  or  higher,  and  we  can  therefore  evaluate 
the  erf  functions  on  the  right-hand  side  of  (A.3)  by  the  stationary  phase  method.  Indeed, 


erf(\/— i?7±) 


0 


fs(u)eiT>±u2du  +  0(6)  , 


o 


where  fs(u)  is  a  specially  chosen  auxiliary  function  such  that  f$  €  C°°[0,  1],  fgk\  1)  =  0  for  all  k  =  0, 1,2, . . ., 
and  fs(u)  =  1  for  u  G  [0, 1  —  <5],  where  <5  >  0  can  be  arbitrarily  small.  Then,  assuming  that  r]±  is  sufficiently 
large,  we  can  apply  the  Erdelyi  lemma  [Fed77]  to  the  last  integral  in  the  previous  formula,  and  by  dropping 
the  0(5)  term  obtain: 


erf  (V~-ir]±) 


fs(u)elri±u~du 


o 


2 r(|) 


-v±- 


rj±_ 

\v±l 


In  doing  so  we  note  that  in  this  particular  case  the  asymptotic  expansion  given  by  the  Erdelyi  lennna  actually 
converges  because  all  of  its  terms  except  for  the  leading  term  happen  to  be  exactly  equal  to  zero.  Altogether, 
if  oj  is  inside  the  interval  [wq  —  B/2,  wq  +  B/2]  and  not  too  close  to  its  endpoints,  then  we  have: 


P(w) 


8  n22y/a  [  |ry_|  |r?+| 

v  > 


1  y/n 
8tt2  y/a 


+*5 

5 


(A.5) 


since  in  this  case  r/_  is  always  negative,  and  //+  is  always  positive. 

Next,  let  u  be  outside  the  interval  [w>o  —  B/2,tu0  +  B/2\  and  again,  not  too  close  to  either  of  its  endpoints. 
Then,  we  can  apply  the  exact  same  reasoning  as  before  and  obtain: 


p  M 


87t2  2y/a 


V- 
\71- 1 


+ 


v+ 
\v+\  J 


=  o, 


=0 


(A.6) 


because  unlike  in  formula  (A.5),  now  both  ?y_  and  77+  have  the  same  sign. 

Finally,  we  need  to  consider  the  case  when  co  is  close  to  one  of  the  endpoints  of  the  interval  [ujq  —  B/2,  wo  + 
B/2],  i.e.,  |  wo  —  w  —  B/2\  or  |w0  -  w  +  B/ 2\  <C  B.  For  example,  if  w  is  close  to  the  left  endpoint  so 
that  |wq  —  to  —  B /2\  is  small,  then  erf (y/—irj_),  see  formulae  (A.3)  and  (A.4),  can  be  evaluated  using  the  first 


AFOSR  Grant  No.  FA9550-07-1-0170.  Final  Performance  Report 


39 


order  Taylor  expansion  for  erf,  whereas  erf(v/—W7+)  should  still  be  computed  with  the  help  of  the  Erdelyi 
lemma,  which  altogether  yields: 


PM 


i  Vn 

87 r2  2,/a 


(wq  —  w  —  B  /  2)  +  1 


(A.7) 


Expression  (A.7)  is  valid  not  everywhere,  but  only  when  the  term  in  the  round  brackets  is  small. 

Hereafter,  we  will  leave  out  the  foregoing  transitional  case,  |wo  —  w  —  B/2\  <C  B  or  |wo  —  w  +  B/2\  <C  B, 
and  evaluate  the  Fourier  transform  P(uj)  only  by  combining  (A. 5)  and  (A. 6): 


P(co) 


ifwe[wo-B/2,wo  +  B/2], 

0,  if  w  ^  [wq  —  B/2,ojq  +  B/ 2]. 


(A.8) 


Qualitatively,  formula  (A.8)  implies  that  in  the  spectrum  of  the  pulse  P(t),  see  (A. 2),  we  have  only  those 
frequencies  u>  that  are  built  in  there  by  design:  wo  —  B/2  <  to  <  loo  +  B /2.  In  reality,  jump  discontinuities 
that  this  pulse  has  at  ±r/2  will  alter  its  spectrum  to  some  degree,  but  we  do  not  take  it  into  account  because 
we  disregard  the  transitional  behavior  (A.7). 

In  the  transformed  space,  each  of  the  frequencies  u>  that  compose  the  spectrum  in  (A.8),  or,  rather,  each 
of  the  waves  elu,t,  propagates  with  its  own  phase  velocity.  Indeed,  the  propagating  pulse  <p(f,r)  that  has 
covered  the  distance  r  from  the  origin  can  be  written  as  the  inverse  Fourier  transform  [BW99,  Ch.  I],  which 
immediately  leads  to  the  introduction  of  the  phase  velocity  uph  =  t'Ph(w)  =  u)/k  in  the  exponent: 


cl>o+B/2  ujq-\-B /2 

(t>(t,z)=  J  P(uj)ei(u,t-kr)duj  =  j  P^y^-r/^^du. 

(jJq  —  B  j  2  OJq — B/2 


(A.9) 


The  function  P(w)  in  formula  (A.9)  is  defined  by  (A.8),  and  the  frequency  w  and  the  wavenumber  k  are 
assumed  to  be  connected  via  the  dispersion  relation  (2.3).  Next,  recalling  that  the  central  carrier  frequency 
is  much  higher  than  the  Langmuir  frequency,  u>q  wpe,  and  also  that  the  radar  pulse  is  narrow  band, 
|w  —  wq I  Cwq,  we  can  use  formula  (2.4)  and  obtain  a  linearized  expression  for  the  phase  velocity: 
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as  well  as  for  its  reciprocal,  which  is  needed  in  (A.9): 
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(A.10) 
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Then,  substituting  the  reciprocal  phase  velocity  into  (A. 9),  we  get: 
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where  A  =  w  -  w<>,  a  =  35  +  and  b  =  t  -  §  (l  +  § . 

Let  us  analyze  the  arguments  of  the  erf  functions  in  (A.ll).  First,  we  can  write: 


t  r  w  B  t  St  def  t' 

aB  =  — | - 4 - =  — I - =  — , 

2  c  ojq  ujq  2  2  2 


(A. 12) 


where  St  is  the  pulse  dilation  that  it  undergoes  when  traveling  over  the  distance  r,  and  t' ,  accordingly,  is 
the  new  pulse  length.  Next,  we  can  linearize  the  group  velocity  (2.5)  the  same  way  as  we  have  linearized  the 
phase  velocity,  see  (2.4)  and  (A.  10),  and,  using  the  first  order  Taylor  formula,  write: 


r  1  cu 

C  2  Cdn 
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VSr(uo)' 


(A.13) 


which  means  that  b  is  the  retarded  moment  of  time  that  corresponds  to  the  center  carrier  group  velocity. 
Then,  by  applying  the  same  stationary  phase  argument  as  we  have  used  for  the  analysis  of  (A. 3),  we  can 
conclude  that  if  b  of  (A.13)  is  inside  the  interval  [— r'/2,  r'/2],  see  (A. 12),  then  the  difference  of  the  two  erf 
functions  on  the  last  line  of  (A.ll)  is  equal  to  2,  while  otherwise  this  difference  is  equal  to  0. 

Finally,  for  the  quantity  l/4a  that  is  in  the  exponent  on  the  last  line  of  (A.ll),  we  can  write: 


1 

4a 


1 


B 


-T  _i_  Ar  P* 
~B 


2t  +  4^^- 

c  UJq  uj  0 


B  def  / 
— — r  =  a  . 
2  t' 


(A. 14) 


Equation  (A.  14)  introduces  the  new  pulse  rate  a1  (slope  of  the  linear  upchirp),  which  is  smaller  than  the 
original  pulse  rate  a.  The  relation  between  the  two  can  be  obtained  as  follows: 


1 


a  = 


2r  (l+2 L^B 
B  V  c  Uq  t 


B[  ru2eB\  rWpeB2  5r  B  def  - 

2?(  1-2c^rj=a--c^4^  =  a-Y^  =  a  +  Sa- 


(A.  15) 


When  deriving  (A.  15),  we  could  use  the  first  order  Taylor  formula  because  for  the  typical  values  of  the 
parameters  involved  (see  the  end  of  Section  2.2  and  the  beginning  of  Section  2.3),  the  second  term  in  the 
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expression  for  a  appears  much  smaller  than  the  first  term.  Indeed,  ~  2.5  -10  12 ,  whereas  for  the 
propagation  distances  r  about  1000/cm  we  have  ~  3.3  •  10~16. 

We  should  also  note  that  even  though  according  to  formulae  (A.  14),  (A.  15)  there  is  a  small  difference 
between  l/4a  and  a,  this  difference  can  be  disregarded  when  computing  the  first  factor  on  the  last  line  of 
(A. 11),  because  taking  it  into  account  would  only  imply  a  minute  change  in  the  amplitude.  Consequently, 

1  _  1 
87t 


and  altogether,  combining  formulae  (A.11)-(A.14)  and  recalling  that  <p(t,r)  =  <f>(t,r)/r,  we  obtain  the 
spherically  symmetric  propagating  pulse  in  the  form: 


<P{t,r)  = 

def 


_Leiu,0(t-r/vphM)XT,(t  _  r/vgr(ujo))eia'(t-r/v^M)2 
4nr 

— A'(f  -  r/vgr(LUo))eiu,o{t-r/v*hM). 


(A.16) 


Appendix  B.  Travel  times  in  the  deterministic  inhomogeneous  ionosphere.  In  this  section, 
the  quantities  Ne  =  Ne(h)  and  wpe  =  uipe(h)  are  assumed  to  have  no  random  component. 


Fig.  B.l.  Schematic  waves’  travel  paths  between  the  antenna  and  the  target  in  the  inhomogeneous  ionosphere. 


Suppose  that  the  antenna  is  positioned  at  x  =  0  and  has  altitude  H  above  the  Earth.  It  sends  a  signal 
toward  the  target  z  on  the  ground,  see  Figure  B.l(a),  so  that  the  look  angle  to  the  target  is  equal  to  90  and 
as  such,  Rz  =  H/  cosf?o-  As  the  medium  is  lossless  (collisionless  dilute  plasma  with  no  Ohm  conductivity, 
see  Section  2.2),  its  electric  permittivity  is  real: 


e(h)  =  1  - 


where  wpe(/i) 


4:TTe2Ne(h) 

me 


(B.l) 


and  we  can  write  Snell’s  law  (for  plane  waves)  in  the  continuous  form  [Gin64,  Sec.  19]  as  follows: 


n(h)  sin  9(h)  =  n(H )  sin  90, 


(B.2) 
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where  the  refraction  index  is  given  by 


n{h)  =  y/eihj  =  \j  1  - 


(B.3) 


The  angle  9{h)  in  formula  (B.2)  is  the  angle  that  the  tangent  to  the  wave  trajectory  at  a  given  altitude  h 
makes  with  the  negative  ordinate  axis,  see  Figure  B.l(a).  Then,  at  every  point  on  the  trajectory  we  can 
write  with  the  help  of  equation  (B.2): 


n(H)  sin  d0 


dx  n(h)s\n9{h)  n(H)  sin d0 

___tan  -  "  n{h)  cos  0(h)  ~~  n(h)  cos  9(h)  ~  ~  ^  {h)  _  n2  {H)  =2  ^  ' 
Hence,  the  actual  trajectory  can  be  obtained  by  integrating  equality  (B.4): 

n{H)  sin  do 


(B.4) 


x(h)  =  /  C^-dh  =  [ 


dh 


H 


n2(h)  —  n2{H)  sin2  d0 


dh. 


(B.5) 


However,  for  the  general  dependence  Ne  =  Ne(h),  the  integration  in  (B.5)  is  not  easy  to  perform.  On  the 
other  hand,  we  can  simplify  formula  (B.4)  by  taking  into  account  that  all  the  ionospheric  corrections  are 
small,  and  employing  the  first  order  Taylor  expansion  with  respect  to  w2e/w2: 


dx 

dh 


1-  sin  d0 


o(H)\ 


sin  dn 


sj  1  -  -  (l  ~  sin^do  \/cos2  d0  + 


^pe  (ff)sin2  e0-ui^(h) 


(B.6) 


-  tan  d0  1  -  - 


1  Wp  e(H) 


2  uj 
\  /  \ 

The  integration  of  equation  (B.5)  yields: 


x  _  1  w_p e(H)  sin"  d0  -  w2e(/i) 
2 


OJ2  COS2  do 


—  tan  do  1  —  - 


1  vL(H)  -  w2(/i) 


OJ2  cos2  d0 


x(h,90)  =  /-tand, 

H  v 

=  (H  —  h )  tan  d0  + 


1  Wpe(#)  -  W2e(/l) 

2 


UJ2  COS2  do 


I  _  *  pe 

1  47re2  tan  d0 


dh 


H 


(B.7) 


2  me  lo2  cos2  do 


J  Ne(h)dh  —  (H  —  h)Ne(H) 


The  shape  of  the  curve  (B.7)  is  determined  by  the  profile  of  the  electron  number  density  Ne  =  Ne(h).  The 
curve  is  also  parameterized  by  the  look  angle  do  at  the  location  of  the  antenna.  With  no  variation  in  the 
electron  number  density,  Ne(h)  =  Ne(H ),  formula  (B.7)  yields  a  straight  line  between  the  antenna  and  the 
target: 


x{h)  =  (H  —  h)  tan  d0. 


(B.8) 


As,  however,  Ne(h)  is  not  constant,  the  ray  that  originates  at  the  antenna  under  the  look  angle  do  will  not, 
generally  speaking,  come  to  the  target,  see  Figure  B.l(a).  To  have  this  ray  come  to  the  target,  we  need  to 
correct  the  look  angle,  see  Figure  B.l(b). 

Let  di  be  the  new  look  angle.  We  substitute  it  into  formula  (B.7)  instead  of  do,  and  require  that 
a;(0,  di)  =  H  tan  do-  In  other  words,  we  require  that  the  wave  trajectory  that  originates  at  the  antenna  under 
the  new  look  angle  9i  terminate  precisely  at  the  target  z  on  the  ground,  see  Figure  B.l(b).  We  have:10 


a;(0,di)  =  H  tandj  + 


1  47re2  NiH)  -  Ne(H)H 

2  me  uj2  cos2  di 


tan  9 1  =  H  tan  do, 


(B.9) 


10Formula  (B.9)  corrects  formula  (D.2)  from  our  previous  publication  on  the  subject  [Tsy09a]  (page  177),  in  which  we  have 
committed  an  arithmetic  error. 
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_  /  rr\ 

where  Ne  is  an  important  characteristic  of  the  ionosphere  —  the  integral  of  its  electron  number  density 
across  the  layer  of  thickness  H: 


JVe' 


u 

Formula  (B.9)  is  an  equation  for  9±.  Let  us  now  introduce  a  new  notation  K  for  brevity: 


K  =  - 

2  TO, 


47re2  iVe(ff)  -  Ne(H)H 
TOe  LO2 


(B.ll) 


as 


and  recast  equation  (B.9) 

K 

Htan6i-\ - —  tan  0i  =  H  tan0o.  (B.12) 

cos2  0 1 

Then,  assuming  that  9\  =  9q  +  5g,  where  59  <C  0o,  we  can  use  the  first  order  Taylor  formula  and  write: 

tan 9\  k,  tan 0O  H - and  — s=s — t--(l  +  2tan0o<50). 

cos 2  0o  cos 2  9 1  cos2  90 

Substituting  these  expressions  into  equation  (B.12),  we  arrive  at 

+  +2t“9”^  +  777)  =  °- 

In  the  previous  equation,  we  keep  only  linear  terms  with  respect  to  59: 

H59  +  K  tan  90  +  x2sm  6°  +  159  =  0, 
cos2  90 


which  yields: 
59  =  - 


K  tan  9 a 


H  +  K  2sin2/»+i 

cos^  6 o 


K  tan  do  1 

ff  1  ,  K  2  sin2  #0  +  1 
1  +  H  ffSgplg 


K tan 0q  /\2sin20o  +  l' 


H 


H  cos2  0o 


(B.13) 


because  according  to  formula  (B.ll),  the  quotient  K/H  is  small,  K/H  ~  w2e/w2  <C  1.  Our  subsequent 
;  shows  that  the  second  term  in  parentheses  on  the  right-hand  side  of  (B.13)  can  be  disregarded. 
s  denote  the  arc  length  of  the  pulse  trajectory.  Then,  using  formula  (B.4)  and  keeping  only  first 
:rms  in  the  Taylor  expansion,  we  can  write: 


ds 

~TT  ~  1/1  + 


2  dx' 
,dh. 


I  ri2(h)  ± 

V  n2(h)  —  n2(H)  sin2  0O  cos0q 


lU2pc(h)  <(g)tan2  0q 

2  U!2 


dh 

T  '  I  •  •  '  / 

Consequently,  the  total  length  of  the  trajectory  which  originates  at  the  angle  9\  and  is 
is  given  by 

1  47re2 


(B.14) 


S  = 


dh 


-  I  H+ - 

cos  0i  \  2  me 


N^H)  -  ule{H)H  2n  \  H  (  K  2 

- - - tan  9 1  = - —  1  +  —  tan  6 1 

/  cos^i  V  H 


LUZ 


is  shown  in  Figure  B.l(b) 

(B.15) 


K  ,, 


0 

For  the  trigonometric  functions  in  formula  (B.15),  we  can  use  (B.13)  and  write: 

-  -  r  ~w~  i  O  ^  /  t  r  /-X  .  9  /I  2^ 


1 


1 


COS  01  COS  00 


n  ,  ±  a  xa\  1  l\  K  tan2  60  f ,  K  2  sin2  0O  + 

(1  +  tan0o<50)  «  - —  1 - - -  I-77 - yv — 

cos  0o  L  H  \  H  cos2  0o 


1 


COS  00 


1  - 


K  tan2  0q 


H 
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and 


t,„1 2  <J,  »  tan2  9„  +  2-^60  *  tan2  «„  -  K  2  ta"2 


cos2  0n 


H  COS2  0n 


1  - 


K  2  sin2  0O  +  1 
H  cos2  0n 


tan2  0o  — 


iF  2  tan2  do 
H  cos2  0o 


On  the  right-hand  side  of  the  previous  equalities,  we  have  dropped  the  terms  of  order  K2/H 2,  because 
K/H  <C  1,  and  hence  the  quantities  quadratic  with  respect  to  K/H  are  negligible.  We  could  have  arrived  at 
the  same  result  if  we  dropped  the  second  term  in  parentheses  on  the  right-hand  side  of  the  last  equality  in 
formula  (B.13)  ahead  of  time.  Substituting  the  previous  two  expressions  back  into  formula  (B.15),  we  have: 


S  = 


H 


cos  0o 


a  K  2 

1  —  —  tan2 

H 


1  +  jj  (  tan2  0O  - 


K  2  tan2  0o 
H  cos2  d0 


H 


COS  00 


K2 

1  +  0|Je 


(B.16) 


Formula  (B.16)  implies  that  up  to  the  negligibly  small  additive  terms  of  order  K2 /H2,  the  length  of  the  pulse 
trajectory  between  the  antenna  x  and  the  target  z ,  see  Figure  B.l(b),  is  equal  to  the  length  of  a  straight 
line  that  connects  the  antenna  and  the  target: 


5: 


H 

COS  00 


=  \z  —  x\  =  Rz. 


(B.17) 


The  travel  time  along  the  trajectory  (B.5)  is  given  by  the  integral 


T_  f  —rlh  -  f  1  ds 

J  dh  J  v(h)dh 


dh. 


(B.18) 


where  v(h)  =  v(u,ujpe(h))  can  be  either  the  group  velocity  (2.31)  or  the  phase  velocity  (2.32),  and  the 
quantity  ||  is  given  by  equation  (B.14).  For  the  case  of  group  velocity,  formula  (B.18)  yields: 


Tgr(x,  z ,  w) 


^pe  V 


1  W2e(0.)  -  W2e(Ff)  ^  2 


costal 


,2 

i+_  p. 


tan  6 i  dh 


o 


1 


H 


c  •  cos  0\ 


l^peW  ,  1  wpeW  -  Wpe(i?) 
2' 


+  7PeV"/tan201  dh, 


(B.19) 


because  we  are  disregarding  all  the  terms  of  higher  order  than  linear  with  respect  to  w2e/w2.  Formula  (B.19), 
along  with  formulae  (B.15)  and  (B.17),  imply  that 


H 


7gr  (o? ,  2 , 02)  ~  + 


_1_  _ 

1 

f  lWpe(/l) 

dh 

1 

c  • 

COS  01 

J  2  W2 

0 

H 

-  -1 — 

1 

1  - 

K  tan2  0( 

C  •  COS  00  C  •  COS  00 


H 


1  47re2  N^H) 


H 

C  ■  COS  00 


1  + 


1  47re2 

2  mew2  H 


c 


1  + 


1  47re2  Nc11 

2  meio2  H 


(B.20) 


where  we  have  again  dropped  the  terms  higher  than  first  order  with  respect  to  K/H  ~  w2e/w2.  Formula 
(B.20)  yields  a  final  expression  for  the  group  travel  time  of  a  high  frequency  narrow  band  radar  pulse  between 
the  antenna  and  the  target  in  an  inhomogeneous  deterministic  ionosphere.  Let  us  also  recall  that  in  our 
linearized  framework  the  length  of  the  pulse  trajectory  is  equal  to  that  of  a  straight  segment  between  the 
antenna  and  the  target,  S  =  Rz,  because  all  the  corrections  are  quadratic,  see  formulae  (B.16)  and  (B.17). 
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Hence,  formula  (B.20)  could  have  been  obtained  merely  as  [cf.  formula  (B.18)] 

R  H 

f  ds  1  f  dh 

^gr (x?Z,Cd)  =  /  =  —  /  y— 

J  Vgr  (h)  cosOo  J  vgr(h) 

0  0 

„  1  1  ( 1  ,  lujUhA  ,,  Rz  (.  ,  1  47re2  7Ve(g)\ 

~c- cos  do  7  (  2  w2  j  c  ^  2mew2  H  )' 


(B.21) 


Formula  (B.21)  can  be  interpreted  in  the  context  of  the  perturbation  theory.  Namely,  the  unperturbed  pulse 
trajectory  is  a  straight  line.  To  get  the  first  order  perturbation  for  the  group  travel  time,  we  integrate  the 
reciprocal  of  the  perturbed  group  velocity  along  the  unperturbed  trajectory. 

As  concerns  the  phase  travel  time,  we  just  need  to  notice  that  the  linearized  group  and  phase  velocities 
differ  only  by  the  sign  of  the  ~  w2e/w2  correction,  see  formulae  (2.31)  and  (2.32).  Hence,  we  have: 


Tph(x,z,u) 


Az  /  1  4tt e2  NjH)\ 

c  y  2  meio2  H  J 


(B.22) 


Appendix  C.  Travel  times  in  the  stochastic  ionosphere.  The  analysis  of  Appendix  B  indicates 
that  we  can  use  the  geometrical  optics  perturbation  method  (see  [RKT89b,  Ch.  I]  for  more  detail)  to  derive 
the  travel  times  of  radar  pulses  in  the  turbulent  ionosphere  as  well.  In  this  case,  the  electron  number  density 
is  a  quasi-homogeneous  random  field  [see  formula  (2.7)]: 

Ne  =  (Ne(h))  +  M(x),  (C.l) 

and  the  dependence  of  fi  on  x  in  formula  (C.l)  means  that  p  depends  on  all  spatial  coordinates,  x  £  M3. 

Next,  with  no  loss  of  generality,  we  can  assume  that  the  propagation  plane  (x,  h)  is  fixed,  see  Figure  B.l, 
and  consider  only  two  independent  spatial  variables.  Moreover,  recall  that  the  random  contribution  to 
the  electron  number  density  is  generally  small  compared  to  the  deterministic  part,  because  the  quantity 
M  =  sj (/x2) / (Ne)  introduced  by  formula  (2.9)  is  small,  see  [Arm05].  Consequently,  we  can  still  linearize 
with  respect  to  the  terms  of  order  Wpe/u;2  as  done  in  Appendix  B  (high  frequency  regime),  and  write  down 
the  equation  for  the  pulse  trajectory  similar  to  equation  (B.6),  but  with  /.t  =  h)  taken  into  account: 

dx_  (  1  47re2  ((Ne(H))  -  ( Ne(h )))  +  (p(0,  H)  -  pQr,  h)) 

dh  0  \  2  TOe  UJ2  COS2  6*0 

Equation  (C.2)  is  a  first  order  ordinary  differential  equation  that  will  be  solved  (approximately)  by  the 
perturbation  method.  First,  we  represent  the  solution  in  the  form  of  a  series: 

x(h)  =x(°\h)+xW(h)  +  ..., 

assuming  that  a:^0^  ~  1  and  ~  w2e/w2,  where  w2e/w2  is  the  small  parameter  for  the  asymptotic  expansion. 
Then,  from  equation  (C.2)  we  immediately  obtain  the  zeroth  order  contribution  to  the  solution: 

h 

x(°\h)  =  J  —  tan90dh  =  (H  —  h)  tan0o,  (C.3) 

H 


(C.2) 


which  is  a  straight  line  between  the  antenna  and  the  target,  see  Figure  B.l  [also  see  formula  (B.8)].  For  the 
first  order  correction,  equation  (C.2)  yields: 


dx^  a  147re2((Ne(H))-(Ne(h)))  +  ^(0}H)-^x^(h),h)) 

-  —  tan  0 o  -  - 


2  me 


UJ2  COS2  Oo 


dh 


(C.4) 
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Note  that  in  order  to  write  equation  (C.4),  we  need  an  additional  assumption  that  the  derivative  of  fi  with 
respect  to  x  be  no  larger  than  order  1.  Integrating  equation  (C.4)  between  H  and  h,  we  obtain: 


c(1)(/i)  = 


1  47re2  tan  0O 

2  me  uj2  cos2  0 


H 


(Ne(h))dh  —  (H  —  h)(Ne(H)) 


H 


(C.5) 


+  f  J  Ai(x(0)  (h),  h)dh  -(H-  h)n( 0,  H ) 


If  /z  =  0,  then  combining  the  zeroth  order  contribution  (C.3)  with  the  first  order  correction  (C.5),  we  arrive 
at  the  previously  obtained  deterministic  solution  (B.7).  Otherwise,  for  the  solution  x(h)  =  x^°\h)  +  x^\h) 
given  by  (C.3),  (C.5)  but  with  the  adjusted  look  angle  9\  we  can  write  similarly  to  (B.9)  and  (B.12): 


:(0, 9i)  =  H  tan  6>i  + 


-  /  rj) 

where  the  quantities  K  of  (B.ll)  and  Ne  ’  of  (B.10)  are  redefined  as  follows: 

/  h 


K 

COS2  01 


tan  0i  =  H  tan0o, 


K  = 


1  47re2 

2  meui2 


(iVe(ff)  -  (Ne(H))H)+i  [  ^x{0\h),h))dh-  n(0,H)H 

^ . - '  J 


deterministic  part 


stochastic  part 


and 


H 


J (Ne(h))dh. 


(C.6) 


=  (K)+rj 


(C.7) 


(C.8) 


Since  sj ( r f)  <  ( K )  or  even  s/W)  «  (K)  [see  the  discussion  right  after  equation  (2.9)],  and  ( K)/H  <C  1, 
we  can  employ  the  same  argument  as  in  Appendix  B  and  show  using  equation  (C.6)  that  for  any  particular 
realization  of  fi  the  look  angle  increment  <50  =  0i  —  0o  will  be  given  by  the  same  expression  (B.13),  but  with 
K  and  N defined  by  (C.7)  and  (C.8),  respectively. 

The  arc  length  along  the  pulse  trajectory  with  the  look  angle  0i  can  be  found  by  integrating  the  differ¬ 
ential  equation  [cf.  formulae  (B.14),  (B.15)]: 

ds  =  1  /  1  4?re2  ((Ne(H))  -  (Ne(h)})  +  -  n(x,h))  2 

dh  cos  0i  \  2  TOe  uj2 


To  integrate  equation  (C.9)  we  can  again  use  the  perturbation  method.  Representing  s(h)  in  the  form 

s(h)=s^(h)+s^\h)  +  ..., 

we  have  the  zeroth  order  term: 

H  -h 


S(°\h)  = 


COSC^I 


(C.10) 


and  the  first  order  term: 


S(1)(M  =  14^e2  tan2  01 
2  me  to2  cos  0i 


H 


(Ne(h))dh  —  (H  —  h)(Ne(H)) 


L  N  h 


H 


+  f  J  n(xt°\h),h)dh-  (H-h)ii(0,H) 


(C.ll) 
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To  obtain  the  full  length  of  the  trajectory  we  substitute  h  =  0  into  formulae  (C.10)  and  (C.ll),  and  also  use 
formula  (C.7),  which  yields: 


S  =  s<°>  (0)  +  (0)  =  ( 1  +  4  tan2  0^  . 

COS  U\  y  ±1  J 


(C.12) 


Expression  (C.12)  formally  coincides  with  (B.15),  only  the  definition  of  K  is  different.  Consequently,  we  can 
repeat  the  argument  from  Appendix  B  and  show  that  for  every  particular  realization  of  p  the  total  length  S 
of  the  pulse  trajectory  will  be  equal  to  that  of  a  straight  line  between  the  antenna  and  the  target,  up  to  the 
terms  of  order  (. K)2/H 2  that  are  negligible,  see  formulae  (B.16)  and  (B.17): 


5  = 


H 

cos  do 


H 


cos  6<o 


=  RZ. 


(C.13) 


Finally,  to  obtain  the  group  travel  time,  we  need  to  integrate  the  following  ordinary  differential  equation: 

dt  1  ds  1  /  1  47t e2  (Ne(h)}  +  n(x,h) 

dh  vgr(x,h)  dh  c  \  2  me  w2 

.  1  (l+-  4?re2  ((^(g))  -  (Ne(h)))  +  (/t(0,  H )  -  fi(x,  h))  tan2 

cos  \  2  me  w2 

1  f,  147re2  (Ne(h))  +  1  4?re2  ((Ne{H))  -  (Ne(h)))  +  (/r(0,  H)  -  /i(z,  h))  2 

c  •  cos  0i  \  2  TOe  w2  2  me  w2  &n  1 

As  according  to  formula  (C.13)  the  shape  of  the  trajectory  will  not  contribute  to  the  travel  time  in  the 
linearized  framework,  the  integration  of  equation  (C.14)  by  means  of  the  perturbation  method  yields: 


0i 


(C.14) 


2gr  (  x ,  z,ui)  — 


H 


C  ■  COS  00 


1  47re2  / 

^  2  ?TigW2  \  H 


H 


+  J  fi(x('°\h),h)dh 


R z_ 
c 


^  1  47re2  (  N( 


\t(h) 


2  mPu!2  \  H 


l  z  N 

+  J  V(s)ds 


(C.15) 


When  deriving  formula  (C.15),  we  took  into  account  that  the  look  angle  increment  is  small,  69  ~  ( K)/H , 
as  we  did  when  deriving  formula  (B.20).  We  also  note  that  the  last  integral  in  formula  (C.15)  shall  be 
interpreted  as  taken  along  the  unperturbed  straight  trajectory. 

The  phase  travel  time,  as  before,  is  obtained  by  changing  the  sign  [cf.  formula  (B.22)]: 


Tph(x,  z,w) 


Rz  f  1  47re2  ( MH) 

c  2  meuj'2  l  H 


1 

R~z 


n(s)ds 


(C.16) 


Appendix  D.  Kirchhoff  integral  and  the  Lorentz  transform.  The  fundamental  solution  of  the 
d’Alembert  operator  is  given  by: 


£(t,z) 


R(t)  6{\z\  —  ct) 
47T  t 


where  R{t)  is  the  Heaviside  function,  and  <S(|z|  —  ct)  is  a  single  layer  of  unit  magnitude  on  the  expanding 
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sphere  of  radius  ct  centered  at  the  origin.  Then,  solution  of  equation  (3.34a)  is  given  by  the  convolution: 

OO 

S(\z  -  z'\  -  c{t  -  t')) 


=  i  /  dt' 


t-t' 


P(t')5(z[  -  vt')S(z'2  -  H)S(z'3)dz' 


1  f  5(\z  —  x(t')\  —  c(t  —  t')) 


47T 


t-t' 
jz  —  x{t')\  d(/i  —  ct) 


P{t')dt' 


(D.l) 


47T  J  ( c\z  —  x(t')\  —  ( Zi  —  vt')v)(t  —  t ') 


P(t')dn 


\x\ 


cP(t') 


47T  c\z  —  x{t')  |  —  {z\  —  vt')l 


fl—Ct 


where  the  new  integration  variable  is  defined  as  /z  =  \z  —  x(t')\  +  ct'  =  \/ ( Zi  —  vt')2  +  ( z^  —  H )2  +  z2  +  ct' . 
To  evaluate  the  final  expression  in  (D.l),  one  needs  to  find  the  value  of  t'  for  which  fi  =  ct.  This  amounts 
to  solving  the  equation 


^ J  (zi  —  vt’)2  +  [z-i  —  H)2  +  z3  +  ct’  =  ct 


(D.2) 


with  respect  to  t' .  Equation  (D.2)  is  a  quadratic  equation,  which  is  very  similar  to  equation  (3.18).  Its 
solution  determines  the  retarded  moment  of  time  when  the  trajectory  of  the  moving  antenna  intersects  with 
the  lower  portion  of  the  characteristic  cone11  that  has  its  vertex  at  ( t.,z ).  The  appropriate  root  of  this 
equation  that  defines  a  retarded  moment  of  time  is 


t'  = 


P2 


t  - 


vz  i 


zl  +  ^-HY  +  zl 


=  _  J_  /(- 1~^)2  .  /  Try,  ..2 _ i  _  \/Cl2  +  (z2  ~H)2+zl 

p  Pc]j  p2  ’  +"3  p  pc 


P  pc 


(D.3) 


where  we  have  used  the  definition  (3.33)  of  the  Lorentz  transform.  Formula  (D.3)  for  t'  is  to  be  substituted 
into  the  numerator  of  the  last  expression  in  (D.l).  As  for  the  denominator  of  this  expression,  using  fj,  =  ct 
we  can  write: 


C\J  ( Zi  —  vt’)2  +  ( Z2  —  H)2  +  z2  —  {z\  —  vt')v  =  c2(t  —  t')  —  ( Z\  —  vt')i 
=  c2t  —  Z\V  —  t'(c2  —  v 2)  =  c2  (t.  —  —  P P 2) 


(D.4) 


=  c 


(2  »;+(«-  Eizzi)  =  ncP) 


where  we  have  used  formula  (D.3).  Altogether,  substituting  (D.3)  and  (D.4)  into  (D.l),  we  see  that  the 
expression  for  ip(t,  z)  reduces  precisely  to  what  is  given  by  formula  (3.35),  as  expected. 

Note  that  equation  (D.2)  that  defines  the  retarded  moment  of  time  t'  is  so  simple  (quadratic)  only  for 
the  case  of  a  uniform  and  straightforward  motion  of  the  antenna.  When  the  motion  is  accelerated  and/or 
non-straightforward,  the  resulting  counterpart  of  equation  (D.2)  may  become  far  more  complicated;  often,  it 
can  only  be  solved  numerically  as  done,  e.g.,  in  [Tsy03].  In  the  theory  of  electromagnetism,  the  corresponding 
solutions  of  the  d’Alembert  equation  can  be  interpreted  as  field  potentials  due  to  moving  charges,  and  in 
this  case  they  are  referred  to  as  the  Lienard-Wiechert  potentials,  see  [LL75,  Chapter  8]. 

11In  the  space  R3,  the  trajectory  is  a  horizontal  straight  line  at  altitude  H ,  see  Figure  B.l.  However,  its  intersection  with  the 
characteristic  cone  needs  to  be  considered  in  the  3  +  1-dimensional  (Minkowski)  space-time,  in  which  the  trajectory  becomes  a 
straight  line  with  the  slope  determined  by  the  velocity. 
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Equation  (D.2)  can  also  be  solved  approximately,  as  done,  e.g.,  in  [CB08].  As  the  square  root  on  its 
left-hand  side  is  the  distance  between  the  observation  point  z  =  (zi,  Z2,  Z3)  €  M3  (location  of  the  target)  and 
the  moving  antenna  at  the  moment  of  time  t' ,  we  can  represent  this  distance  using  the  law  of  cosines  and 
then  derive  an  approximate  expression  with  the  help  of  the  Taylor’s  formula.  Let  7  be  the  angle  between 
the  velocity  v  and  the  direction  from  the  antenna  x{t')  to  the  target  z ,  see  Figure  B.l;  also  recall  that 
r  =  \Jz\  +  (Z2  -  H )2  +  z2.  Then, 

yj (zi  —  vt')2  +  (z2  —  H)2  +  z2  =  \/r2  +  (vt’)2  —  2 rvt’  cos  7 


/  2 rvt'  cos  7  /  vt'\ 

=  r\  1 - -  +1  —  I 


r  J 


'■  r  —  vt'  cos  7. 


In  this  case,  equation  (D.2)  becomes: 


r  —  vt'  cos  7  +  ct'  =  ct 


so  that  its  solution  is 

,  t-  ^ 

t'  =  _ £ _ 

1  -  cos  7 

Substituting  this  expression  into  formula  (D.l),  we  arrive  at  the  solution  (3.40). 
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